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Abstract 

We address the problem of the response to a seismic wave of an urban site consisting of N 
non-identical, non-equispaced blocks overlying a soft layer underlain by a hard substratum. The 
results of a theoretical analysis, appealing to a space- frequency mode-matching (MM) technique, 
are compared to those obtained by a space-time finite element (FE) technique. The two methods 
are shown to give rise to the same prediction of the seismic response for = 1 and N — 2 blocks. 
The mechanism of the interaction between blocks and the ground, as well as that of the mutual 
interaction between blocks, are studied. It is shown that the presence of a small number of blocks 
modifies the seismic disturbance in a manner which evokes qualitatively, but not quantitatively, 
what was observed during the 1985 Michoacan earthquake in Mexico City. Disturbances at a 
much greater level, induced by a large number of blocks (in fact, a periodic set) are studied in 
the companion paper. 
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1 Introduction 



The Michoacan earthquake that struck Mexico City in 1985 presented some particular charac- 
teristics which have since been encountered at various other locations |4 H I42L [77| Ki4 | I25j . but at a 
lower level of intensity. Other than the fact that the response in downtown Mexico varied consid- 
erably in a spatial sense ,15 , was quite intense and of very long duration at certain locations (as 
much as ~3min "^S"), and often took the form of a quasi-monochromatic signal with beatings |36", a 
remarkable feature of this earthquake (studied in 1^ ^ .21, 22_, ) was that such strong motion could 
be caused by a seismic source so far from the city (the epicenter was located in the subduction zone 
off the Pacific coast, approximately 350km from Mexico City). It is now recog nized EJ that the 
characteristics of the abnormal response recorded in downtown Mexico were partially present in the 
waves entering into the city (notably 60km from the city as recorded by the authors of jE]) after 
having accomplished their voyage from the source, this being thought to be due to the excitation 
of Love and generalized-Rayleigh modes by the irregularities of the crust El IS] ) • 

In the present investigation (as well as in the companion paper), we focus on the influence of 
the presence of the built features of the urban site as a complementary explanation of the abnormal 
response: the so-called city-site effect. A building or a group of buildings over a hard half-space, 
solicited by a plane incident SH wave, has been shown to modify the seismic waves on the ground 
near the building |46| I33j . the modification being larger when more buildings are taken into account 
because of multiple-interaction: i.e., the so-called structure- soil- structure interaction. For models of 
the geophysical structure involving only a hard half-space, the stress-free base block mode appears 
to be the main cause of the modification j33j . 

The studies that deal with a geophysical structure involving, in addition, a soft-layer overlying 
the hard-half space, have been mainly concerned either with an infinite set of periodically-arranged 
[H 1^ [3] or randomly- arranged [321 llOj buildings on, or partially imbedded in, the ground. In 
[2], the authors suggest that the large duration and amplitude are strongly linked to resonant 
phenomena of the soft-layer associated with waves whose structure is close to that of Love waves. 
The solicitation being of the form of a plane incident wave, such modes cannot be excited in the 
absence of buildings |U . 

In |3H I26j . it was shown that the modes of a soft layer/hard half space can be excited when 
the interface between the subtratum and the layer present some irregularities. These effects were 
qualified as "vertical and lateral interferences" in a previous numerical study T. The question of 
the excitation of modes, via surface irregularities constituted by the set of buildings on the ground, 
was subsequently addressed in jTH] . In |3H1 it was found that the excitation of vibration modes 
associated with a periodically-modulated surface impedance, modeling a periodic distribution of 
blocks emerging from a flat ground, can lead to enhanced durations and amplifications of the 
cumulative displacement and velocity as compared to what is found for a flat stress-free or constant 
surface impedance surface. The authors of 48_ show that these modes manifest themselves by 
amplified evanescent waves in the substratum. 

The contributions |^ 1^ employ homogeneized models of a periodic city, but the fact that these 
models are restricted to low frequencies may explain why they do not account for the amplifications 
obtained in [TO] [IS] . In a host of other numerical studies [Sl[ll[Hlll[Sl[inilISllIl,the presence of 
buildings is found either to hardly modify, or to de-amplify, the seismic disturbance, in contradiction 
with what is shown in [431 1471 117j . 

In [SH], the spatial variability of damage to structures on the ground was attributed to the 
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variability of the resonance frequencies of the buildings and of the soil structure beneath each 
building, with the implication that the most dangerous situation is when the natural frequency of the 
building (often treated as a one degree (or several degrees) of freedom oscillator [23 ESI EB El EII ) is 
coincident with that (obtained by a ID analysis) of the substructure below the base of the building 
(a well-known paradigm in the civil engineering community known as the double resonance). 

Another point of view is to consider the building as a seismic source, either when it is solicited 
artificially by a vibrator located on its roof |28| \50\ , or when it re-emits vibrations received from the 
incident seismic disturbance (or other form of solicitation such as that coming from an underground 
nuclear explosion |44| [T2j ) . It is not unreasonable to think that the presence of one or more buildings 
on the ground enables the excitation of the (Love, Rayleigh) modes of the underground system. 
This is known to be possible when a flat stress-free surface overlying a soft layer in welded contact 
with a hard substratum is solicited by a source located in the layer or substratum |2J and should 
therefore also occur when the source (i.e., the building) is on the free surface. 

The present work originated in the observation that no satisfactory theoretical explanation 
has been given until now of the influence of buildings on anomalous seismic response in urban 
environments with soft layers, or large basins, overlying a hard substratum. The principal reason 
for this knowledge gap probably lies in the complexity of the sites examined in previous studies 
and in the complexity of the phenomena. Thus, it appears to be opportune to develop a theoretical 
model which is as complete and as simple as possible, on an idealized, although rather representative 
urban site, in order to address the following questions: 

(i) how should one account for the principal features of the seismic response in the cases of a 
relatively small, and then, large number of blocks? 

(ii) what are the modes of the global structures (i.e. the superstructure plus the geophysical 
structure) and what are the mechanisms of their excitation and interaction? 

(iii) what are the repercussions of resonant phenomena on the seismic response? 

(iv) what are the differences in seismic response between configurations with a small and a large 
number of blocks? 

The investigation herein focuses on the seismic response of one and two blocks (the case of a 
periodic set of blocks is considered in the companion parper) in welded contact with a soft layer 
overlying a hard half-space. The modal analysis of the whole configuration, backed up by extensive 
numerical computations, shows that: i) the presence of one block induces the excitation of two 
types of modes, the first whose structure is close to that of a mode of the geophysical structure, 
and the second whose structure is close to that of a mode of the superstructure (i.e., the set of 
blocks, each of which is formed of one or several buildings), ii) the presence of more than one block 
gives rise to coupled modes resulting from a combination of the two other types of modes. 

We uncover the mechanism of the (so-called soil-structure) interaction between the superstruc- 
ture and the geophysical substructure. Despite the fact that differences are noticed between the 
computed displacements for a configuration with, and in the absence of, buildings, mainly consist- 
ing in a longer duration and a larger displacement in the building than at the same location in the 
absence of the building, and in a modification, due to the presence of the block(s) of the structure 
of the waves traveling in the layer, no very pronounced effects are apparent in the case of only a 
few (one or two) buildings. This could mean that the city-site effect is important only when a large 
number of buildings is accounted for. This case is investigated theoretically and numerically in the 
companion paper for a periodic arrangement of identical blocks. 
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Figure 1: View of a modern city with the underground. 



2 Description of the configuration 

We focus on a portion of a modern city consisting of a set of blocks (see e.g., fig. ^in which 
it will be noted that the blocks (i.e., buildings or groups of buildings) are not generally identical, 
nor arranged in periodic manner). The city has 2D geometry, with X3 the ignorable coordinate 
of a 0x1X2X3 cartesian coordinate system (see Fig. EJ. The buildings are assumed to be in 
welded contact, across the flat ground surface, with the substructure. The latter is composed of a 
horizontal soft layer underlain by a hard half space (see fig. Each block is characterized by two 
constants, its height bj and width wj, and all blocks have the same rectangular geometry (but not 
the same sizes) and composition. Let dj be the xi coordinate of the center of the base segment of 
the j— th block. The distance between the blocks j and i is denoted by dji = \dj — di\ and is not 
necessarily constant between successive pairs of blocks. 

For the purpose of analysis, each block is homogenized (this does not mean that the set of blocks 
is reduced to a single horizontal, homogeneous layer, as in [3] ) , so that the final aspect of the city 
is as in Fig. 13 Let B S Z denote the set of indices by which the blocks are identified (e.g., for three 
blocks: {1,2,3} or { — 1,0,1}). The cardinal of B is designated by (i.e., denotes the number 
of blocks in the configuration, and this number will either be finite (in the following analysis) or 
infinite (as in the companion paper). 

Fj is the stress- free surface composed of a ground portion Tg, assumed to be flat and horizontal, 
and a portion Tag, constituting the reunion of the above-ground-level boundaries Tag ; j £ M of 
the blocks. The ground Tq is flat and horizontal, and is the reunion of F^ and the base segments 
F;^^ ; j G B joining the blocks to the underground. 

The medium in contact with, and above, Fj is air, assumed to be the vacumn (which is why Tf 
is stress- free) . The medium in contact with, and below To is the mechanically-soft layer occupying 
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Figure 2: View of the 2D city (only two of the blocks are represented). 




Figure 3: View of the 2D city with homogenized blocks (only two of the blocks are represented). 



Figure 4: Sagittal {xi0x2) plane view of the 2D city with homogenized blocks (only two of the blocks 
are represented) solicited by a cylindrical wave radiated by a line source located at (xf, j;!) S ^o- 
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Figure 5: Sagittal {xi0x2) plane view of the 2D city with homogenized blocks (only two of the 
blocks are represented) solicited by a plane wave with incident angle 9^ . 



the domain Oi, which is laterally-infinite and of thickness h, and whose lower boundary is F/j, also 
assumed to be flat and horizontal. The soft material in the layer is in welded contact across F/j 
with the mechanically-hard material in the semi- infinite domain (substratum) ^Iq. 

The domain of the j-th block is denoted by and the reunion of all the fig > J G B is denoted 
by 0,2- The material in each block is in welded contact with the material in the soft layer across 
the base segments F^^ ; j G B. 

The origin O of the cartesian coordinate system is on the ground, X2 increases with depth and 
X3 is perpendicular to the (sagittal) plane of the figs. I4I5I With ij the unit vector along the positive 
Xj axis, we note that the unit vectors normal to Tq and F/j are —12. 

The media filling Qq, Jli and IJjgB^2 respectively and the latter are 

assumed to be initially stress-free, linear, isotropic and homogeneous (thus, each block, which is 
generally inhomogeneous, is assumed to be homogenized in our analysis). We assume that is 
non-dissipative whereas and are dissipative, described by a constant quality factor in 
the frequency range of excitation. 

The seismic disturbance is delivered to the site in the form of a shear-horizontal (SH) cylindrical 
wave (radiated by a line source parallel to the X3 axis and located in J7o; see fig. ^ or a plane 
wave (with incident angle 0*; see fig. EJ, propagating initially in JIq (this meaning, that in the 
absence of the layer, the city, and the air, the total field is precisely that associated with this 
cylindrical or plane wave). The SH nature of the incident wave (indicated by the superscript i 
in the following) means that the motion associated with it is strictly transverse (i.e., in the ^3 
direction and independent of the 2:3 coordinate). Both the SH polarization and the invariance of 
the incident wave with respect to X3 are communicated to the fields that are generated at the site in 
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response to the incident wave. Thus, our analysis deals with the propagation of 2D SH waves (i.e., 
waves that depend exclusively on the two cartesian coordinates {xi,X2) and that are associated 
with motion in the direction only). 

We shall be concerned with a description of the elastodynamic wavefield on the free surface 
(i.e., on Ff) resulting from the cylindrical or plane seismic wave sollicitation of the site. 

3 Governing equations 

3.1 Space-time framework wave equations 

In a generally-inhomogeneous, isotropic elastic or viscoelastic medium M occupying M^, the 
space-time framework wave equation for SH waves is: 



wherein u is the displacement component in the is direction, / the component of applied force 
density in the 13 direction, fi the Lame descriptor of rigidity, p the mass density, t the time variable, 
uj the angular frequency, 9" the n— th partial derivative with respect to t, and x = (xi,X2). Since 
our configuration involves three homogeneous media, and the applied force is assumed to be non 
vanishing only in Qq, we have 



wherein superscripts m designate the medium (0 for M^, etc.), 6fn0 = 1 for m = 0, 5mo = for 
m / 0, and c™ is the generally-complex velocity of shear body waves in M"^, related to the density 
and rigidity by 



it being understood that p™^, fi'^{uj) ; m = 0,1,2 are constants with respect to x. In addition, 
the densities are positive real and we assume that substratum is a dissipation-free solid so that the 
rigidity therein is a positive real constant with respect to lo, i.e., ix^{uj) = ijP > 0. 

3.2 Space- frequency framework wave equations 

The space-frequency framework versions of the wave equations are obtained by expanding the 
force density and displacement in Fourier integrals: 



V • (/x(x,u;)Vu(x,0) - p(x)a|u(x,0 = -p(x)/(x,0 , 



(1) 



(c'"(u;))2 V • V«'"(x, t) - a2u™(x, t) = -/(x, t)dmO ; X G ; m = 0, 1, 2 , 



(2) 





so as to give rise to the Helmholtz equations 



V • Vu"(x,u;) + (/c™(a;))2^i-(x,a;) = -/(x,a;)<5^o 



Vx G a 



"m 



m = 0, 1 , 



(5) 



wherein 




(6) 
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is the generally-complex wavenumber in M™. Actually, due to the assumptions made in sects. [21 
andO 

fcV):=^=-y|, (7) 

(i.e., is a positive real quantity which depends linearly on w ). 

As mentioned above, we shall be concerned with cylindrical or plane wave excitation of the city. 
Plane waves correspond to / = and cylindrical waves to / 7^ 0. 

The incident field is chosen to take the form of a pseudo Ricker-type pulse in the time domain. 

3.3 Space-frequency framework expression of the driving force for cylindrical 
wave excitation 

The space-frequency framework expression of the driving force density for a cylindrical wave 
radiated from a line source located at x** := {x\,X2) G f^o is 

/(x,a;) = SH5(x-x^) , (8) 

wherein S{oj) is the spectrum of the incident pulse and is chosen to be a time derivative of a Ricker 
pulse and 5{.) is the Dirac distribution. The amplitude spectrum S{uj) is given by 

^M = — 7^i^«-pK--i^) ' (9) 



to which corresponds the temporal variation (Fourier inverse of S{u))): 

S{t) = -247rQ^ [-3(i, -t) + 2a^{ts - tf] exp [-a^{ts - tf] , (10) 

wherein a = vr/tp, tp is the characteristic period of the pulse, and tg the time at which the pulse 
attains its maximal value. In the remainder of this paper, we shall take tg = tp = 2 sec. 
The (incident) wave associated with this driving force is 

u\^,u;)= [ G°(||x-y||,.;)/(y,u;)dt^(y) , (11) 

wherein y := (yi, 2/2) is an integration point in the sagittal plane, dw{y) the differential area element 
at point y and G°(A;°||x - y||) the 2D free -space Green's function which satisfies: 

A+(fcO)']G°(||x-y||,a;) = -5(x-y) ; Vx G M^^ (12) 

(with 6{ . ) the Dirac delta distribution) and the (outgoing wave) radiation condition: 

G'^dlx — y||, w) ~ outgoing waves ; ||x — y|| — > 00 . (13) 
The free-space Green's function is given by |37] : 

G°(||x - y||, ..) = ^Hl,'\k'\\^ - y||) = ^ exp{i[fci(xi - yi) + ^01x2 - 2/2!]}^ , (14) 
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with H^\.) the Hankel function of the first kind and order 0, and 



A:^ = ^{kif - {kif ; > , > for w > . (15) 

Introducing (fTl|) and Q into (|TT|l results in 



ip(i)^;.o 

4 ' 

which is the space-frequency expression of a cyhndrical wave. 



u'{^,uj) = S{uj)GW\j^-^'\\,uj) = S{uj)-Hl,^\k^\\^-Ji.'\\) , (16) 



3.4 Material constants in a dissipative medium 

A word is now in order concerning the dissipative nature of the layer and blocks. In seismolog- 
ical applications involving viscoelastic media, the quality factor is usually considered to be either 
constant or a weakly-varying function of frequency jEl in the bandwith of the source. We shall 
therefore assume that Q^{lo) = , with constants, j = 1,2. It can be shown [^01 that this 
implies 

/ • \ — arctan ( ~ ] 

M'M = f.;.,(^)' *°'^J = 1.2, (17) 

\OJref J 

wherein: w^e/ is a reference angular frequency, chosen herein to be equal to 9 x 10~^Hz. Hence 

^^^) = <eA^\ ;j = i,2, (18) 



re/ 



with c^gj := Y Note should be taken of the fact that even though , j = 1,2 are non- 

dispersive (i.e., do not depend on lS) under the present assumption, the phase velocities e' ; j = 1, 2 
are dispersive. 



3.5 Boundary and radiation conditions in the space-frequency framework 

The translation of the stress-free (i.e., vanishing traction) nature of = Fg |J Fa^, with F^^ := 

^i(w)5„ui(x,w) = ; xEFg, (19) 

/x2(a;)9„u2(^)(x, c^) = ; X G F^^ , j G B (20) 

wherein n denotes the generic unit vector normal to a boundary and (?„ designates the operator 
9„ = n • V. 

That and are in welded contact across F^s := UjeB-'^L translated by the fact that 
the displacement and traction are continuous across F^^: 

M^(x,cj) - m2(^')(x,cj) = ; X G F^^ , j G B (21) 
li^{uj)dnU^{^M - ;u2(cj)5„u2(^')(x,u;) = ; X G F^^^ , j G B. (22) 
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That Afi and M° are in welded contact across F/j is translated by the fact that the displacement 
and traction are continuous across this interface: 

u^(x, w) — u^{'x,u)) ; X G F/i , (23) 

fi\Lo)dnu\^,uj) - fi^{Lj)dnu\^,Lo) ; X G T;,. (24) 

The uniqueness of the solution to the forward-scattering problem is assured by the radiation con- 
dition in the substratum: 

u^{x,uj) — n'(x,Lj) ~ outgoing waves ; ||x|| oo, X2> h . (25) 

3.6 Boundary and radiation conditions in the space-time framework 

Since our finite element method jl8l 1191 117j for solving the wave equation in a heterogeneous 
medium M (in our case, involving three homogeneous components, , and M^) relies on the 
assumption that M be a continuum, it does not appeal to any boundary conditions except on Tj 
where the vanishing traction condition is invoked (fictitious domain method). Furthermore, since 
the essentially unbounded nature of the geometry of the city cannot be implemented numerically, 
we take this geometry to be finite and surround it (except on the Tf portion) by a perfectly- 
matched layer (PML) which enables closure of the computational domain without generating 
unphysical reflected waves (from the PML layer). In a sense, this replaces the radiation condition of 
the unbounded domain. The stress-free boundary condition on Fj is modeled with the help of the 
fictitious domain method jSI , which allows us to account for the diffraction of waves by a boundary 
of complicated geometry, not necessarily matching the volumic mesh. 



3.7 Statement of the boundary-value (forward - scattering) problem in the 
space-time framework 

The problem is to determine the time record of the displacement fields n^(x, t) on Tg and 
^/20)(x,t) on Tig,jeM. 

3.8 Recovery of the space-frequency displacements from the space-time dis- 
placements 

The spectra of the displacements are obtained from the time records of the displacements by 



Fourier inversion, i.e.. 



n^ (x, ^) = ^ / ^^' Cx, t)e''^'dt ; j = 1, 2 . (26) 



4 Field representations in the space- frequency framework for < 

oo 

4.1 Field in VIq 

It is useful to consider the boundary dO,Q of to be composed of plus a semi-circle of 
infinite radius Fqo joining F/^ at x = {—oo, h) and x = (-co, h). The unit vector n normal to dVt^ 
is taken to be outward with respect to VLq, so that it is equal to —\2 on T^. 
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We seek the field representation in Qq. Applying Green's second identity to and in Qq 
and making use of the radiation condition at infinity relative to these two functions, gives 

?t:f7o(x)n°(x,a;) =u'{x,u;) + 

[GO(||x-y|U)a„txO(y,a.)-^/0(y,w)5„GO(||x-y||,u;)]d7(y) , (27) 



- h 

wherein d^^y) is the infinitesimal arc length along and 



1 ; X G 

nno{^) = { ;yGR''\Tlo . (28) 

1/2 ; y G Th 

Introducing the cartesian coordinate integral representation of the Green's function 1)141) into the 
boundary integral representation of the field ()27j) . while paying attention to the absolute values, 
leads to the following result: 

/°° dk^ 
S°(A;i,cj)exp{i [kixi + k^{x2 - h)]} ; x G J^o , (29) 
-oo "-2 

wherein: 



B^{ki,uj) = —J {dy^u°{yi,h,uj)+ik^u'^{yi,h,uj)}exp{-ikiyi)dyi , (30) 

At this point, we must distinguish between plane wave excitation (briefly alluded-to in a subsequent 
section) and cylindrical wave excitation (to which all the following numerical results apply). 
In the case of plane-wave excitation we can write 

A°-{ki,uj)exp{i[kiXi-k^X2]}^ ; x G , (31) 

-oo ^2 

wherein 

A°-{ki,Lo) = S{uj)kl6{ki - k\) , (32) 
with k\ = k^sinO^, and 0* the angle of incidence, so that 



u 



(x, Lo) = S{lo) exp {i [k{xi - k\x2] } ; x G , (33) 



wherein /cg = k^cosO"^. 



In the case of cylindrical wave excitation, we have, on account of (|14|) and ()16|) : 

' f°° dk^ 

/ A^+{ki,uj)exp{i[kiXi + k^X2]} ^ ; x G 0+ 

J -oo ^^2 

f°° dki 
I A^~{ki,uj)eyi-p{\[kixi-klx2]] -r-Q ; x G 

J -oo "-2 



n*(x, w) = < 



(34) 



wherein 

0+ = {Vxi G M ; X2>xl] , (35) 
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Q- = {Vxi G R ; /i< X2 < x|} , (36) 

A°+{ki,uj) = S{uj)^exp{-i[kixl + k^x'2]} , (37) 

A^-{ki,uj) = S{uj)^exp{-i\kixl-k^x'2\] . (38) 
Since we shall henceforth be interested only in the field in the subdomain of Oq, we can write 

/°° dki 
A^~{ki,Lv)exp {i [kixi - k2X2] } -70"+ 
-00 ^2 

B^{ki,uj) exp {i [A:ixi + A:0(x2 - /i)] } ; x e , (39) 

with the understanding that: i) S{uj) is known a priori and given by its expression in @, ii) 
A^^{ki,ui) are known a priori, iii) i?*^(A:i, w) is an unknown function, iv) A^^{ki,u!) and B^{ki,u)) 
have units of (length) since (and, in general, all displacements), have units of length, v) is 
expressed as a sum of incoming and outgoing plane (bulk and evanescent) waves. 

4.2 Field in Qi 

By proceeding in the same manner as previously we find 

/°° dki 
A^{ki^uj) exp |i \kixi — k\x2\ } -r^^ — h 
-00 "-2 

/°° dk^ 
5^(A;i,a;)exp{i [kiXi + k\x2]] ; x G J7o , (40) 
-00 1^2 

with the understanding that: i) now both A^{ki,uj) and B'^{ki , u}) are unknown functions, ii) v} is 
expressed as a sum of incoming and outgoing plane (bulk and evanescent) waves. 

4.3 Field in the j-th block 

The task is here to obtain a suitable representation of the field in the generic block j (of height 
hj and width Wj) occupying the domain The boundary of this domain is = Tagljr^^. It 
should be recalled that the field satisfies a Neumann boundary condition on the emerged boundary 
Tag of the block. No boundary condition of the Neumann or Dirichlet type is available on the 
segment F;^^ so that, strictly speaking, we are not searching for a modal representation of the field 
in the block domain, but rather for a quasi-modal representation, the latter satisfying a priori the 
boundary condition on Tag-, but no particular boundary condition on F;^^. 

Let 0^x\x2X'^'^ be the (local) cartesian coordinate system attached to ^2 such that the origin 

is located on, and at the center of, the segment F^^. We note that 

xi = dj + x{ , X2 = X2 ; Vj G B , (41) 
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wherein it should be recalled that dj is the xi coordinate of the center of the base segment of the 
j—th block. 

We apply the separation of variables technique and the boundary conditions on SJlg to obtain 



,2(i), 



x,a;) = 5]B^(^-)(a;)cos 



m=0 



cos 



x£n{ , yj £M , (42) 



wherein 



and Bm''^ (uj) has units of length. On account of (|4H) we finally obtain 



mvr 



5 k' 



2(i) 
2m 



{kr-{& ; 3^teM>0 , Q{kZ')>0 fora;>0, (43) 



Xi) 



,2(i) 



m=0 



(x,^)=^i?^0)(a;)cos klliUx,-d,+ 



Wi 



COS 



X G , Vi G B , (44) 



it being understood that the B2(j') := {Bm''\io) ; m G Z} , j G B are all unknown vectors. 



5 Determination of the various unknown coefficients by applica- 
tion of boundary and continuity conditions on Vg and Th for the 
case N < oo 

5.1 Application of the boundary and continuity conditions concerning the trac- 
tion on 

From ((201) and JSIl) we obtain 

/•oo 



di+■w/2 
jgB Jdj-w/2 



a^.2M^^^'^(xi,0,a;)exp(-iKiXi)(ixi = ; VETi G M . (45) 



Introducing the appropriate field representations therein and making use of the orthogonality con- 
dition 



/oo 
exp[-i(/ci - Ki)xi]dxi = 2-K5{ki - Ki) ; VA;i , /Ci G M , 
-oo 

gives rise to 



(46) 



27ri ^ 



oo 2 7,2(i) 



m=0 



wherein 



Im^{ki,uj) := / exp{±ikiWjri)cos{klll^Wjr])dr] 
Jo 



) ; Vfci G M , (47) 



(48) 
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5.2 Application of the continuity conditions concerning the displacement on 
From ((221) obtain 



i^(a;i,0,u;)cos klll\xi - di + wi/2) 



dxi — 

u^^^\xi,Q,uj)cos\kll\xi- di + wi/2)\dxi = {) ; V/ G B . (49) 



Introducing the appropriate field representations therein, and making use of the orthogonahty 
condition 



di- 



cos 



^Vrrli^i- di + wi/'^) cos k1^[' {xi - di + Wi /2) 



rl - 

(1X\ — Omn 



; ym , n = 0,1,2,.... , (50) 



wherein 5mn is the Kronecker symbol and Cm the Neumann symbol (e^ = 1 for m = 0, = 2 for 
m > 0), gives rise to 



[A\k,,u:)+B\k,,u;)] {k,,u.)e'^^^'^--^/^)^-^ 



; Vm = 0,1,2,.... , l€M . (51) 



5.3 Application of the continuity conditions concerning the traction on Th 

From (|2,''{|) we obtain 

/OO f'OO 
dx2U^{xi,h,uj) exp{—iKiXi)dxi — fi^ / dx2U^{xi,h,uj) exp{—iKiXi)dxi = ; VETi G M . 
-OO J — OO 

(52) 

Introducing the appropriate field representations therein, and making use of the orthogonality 
relation (|l6|) . gives rise to 



-/iM"-(fci,w)e-i'^"^ + /i°5°(A;i,Lj) + ;uM^(A;i,w)e-^''2^-/S^(fci,a;)ei'^^ft = ; Vfci G M . (53) 

5.4 Application of the continuity conditions concerning the displacement on Th 

From H24p we obtain 

/OO fOO 
u^{xi,h,u>) ex-p{—iKiXi)dxi — / u^{xi,h,u>) exp{—iKiXi)dxi = ; Vi^i G M . (54) 
-OO J —OO 

Introducing the appropriate field representations therein, and making use of the orthogonality 
relation (|l6|) . gives rise to 



/V2 



1,0 
"-2 



^2 



^2 



(55) 
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5.5 Determination of the various unknowns 

5.5.1 Elimination of Bm{uj) to obtain an integral equation for cj) 

After a series of substitutions, the following integral equation for B^{ki,uj) is obtained (wherein 
Ki and K2 play the same roles, and are related to each other in the same manner, as ki and fcg 
respectively) : 

/oo 
D{ki,Ki,uj)B^{Ki,uj)dKi=F{ki,uj) ; VA;i G R , (56) 
-00 

wherein 

C{ki,uj) = cos{klh) - sin{klh) , (57) 



D{kuKi,io) 
i 

2^ 



cos{Klh) - i^^ s\n{Klh) 



^gi(Ki-fci 



2 

00 

m=0 



2, 2{j) 



; V/ci, Ki G M , (58) 



and 



F{ki,Lo) = A^~{ki,uj)e-'''> 



cos{klh) + i Q sin(/i;2/i) 



+ 



cos(E:i/i)+i^sin(Ki/i) 



00 



)(dj-u,j/2) 



m=0 



tan(A:2'(;^)6,)4;?'^'(^i,^)^^'^+(i^i,^)rf^i 



; V/ci G M , (59) 



Eg . (1561) is a Fredholm integral equation of the second kind for the unknown function {B^{ki,uj) ; ki G 
M}. 

By means of the changes of variables ki = k^ai, Ki = k^Si (note that ai and are dimen- 
sionless), we can cast (|56|) into the form 



/oo 
E{ai,Si,uj)B^{Si,uj)dSi = F{ai,uj) ; Vai G M , 
-00 



(60) 



wherein C((Ti,u;), B^{ai,u}) and F{ai,uj) are, by definition, C{ki,uj), B^{ki,Lo) and F{ki,uj) in 
which A;i is replaced by k^ai, D{(Ti, Si,lo) is D{ki, Ki,uj) in which fci, i^i are replaced by k^ai, 
K^Si respectively, and 

E{ai,Si,(jj) ■.= k^D{ai,Si,uj) . (61) 
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Adding and subtracting the same term on the left side of the integral equation gives 

/oo 
E{ai,Si,Lo)[l- S{Si- ai)]B^{Si,io)dSi = F{ai,uj) ; Vai G M , (62) 
-QO 

from which we obtain 

B (0-1,0;) = — r — r ; V(7i e M . (63) 

C[ai,u}) - E{ai,ai,u)) 

An iterative approach for solving this integral equation consists in computing successively: 



S0(0)( 



C{ai,Lo) - E{ai,ai,Lo) 



; VcTi G R , 



(64) 



50W(ai,u;) = i?(°) 



J^^E{ki,Si,iv) [1 - SjSi - ai)] i?°W(gi,o;)d5i 
C((Ti,w) - E{ai,ai,u) 



■ V(7i e M , (65) 



and so on. 



5.5.2 Elimination of B^{ki,uj) to obtain a linear system of equations for fim ('^) 

The procedure is again to make a series of substitutions which now lead to the linear system of 
equations for bI^^\u), V/ G B, Vn G N : 

C,(o;)52(0(^) = + ^Y1 D^nil{^)B^l'\^) ; VI G B ; Vn G N , (66) 

j& m€N 



wherein 



Cn{oj) = cot(A;i„6) ; Vn G N 



sm(A;2„ 0/) -^-oo 



4^^+(fci,a;)e^^i('^'-""/^) 
cos{klh) — i^^ow sin{klh) 



dki 
To" 



and 

^ ^ ifeSV^^^j^n Sin(fc2^^bj) ^^ 
27rsin(4')60/^° 

/oo cos(klh) — i^Vrr sin(A;9/i) 

I«+(fci,a;)/^)-(^i,.) \l' 
-00 



cos I 

cos{klh) - ij^sin(A;|/i) 



(67) 



; V/ G B ; Vn G N , (68) 



ik,{(d,-d,)-^)dki 
VI, j G B ; Vn, m G N . (69) 
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Eq. can be written as: 



B^^'->{uj) = jjj- 

Cn{u:)-D^Jirl{u:) 

; V/ e B ; Vn e N , (70) 
An iterative procedure for solving this linear set of equations is as follows: 

Bf^ (c.)) = ^ , f^^^l , , ; yieM;yneN, (71) 



{bi^^Hu^: 



C„(a;)-Di'i)(c^) 

; p = 1,2, ... ; V/ € B ; Vn € N . (72) 



6 Modal analysis 

6.1 The emergence of the natural modes of the configuration from the iterative 
solution of the integral equation for B^(ki,uj) 

Eqs. ()64() . ()65() . etc. show that the n-th order iterative approximation of the solution to the 
integral equation (|5H|) is of the form 

C{ai,uj) - E{ai,ai,LO) V{ai,LO) 

wherein 

Ar(o)(^7i,c^) = F(ai,c^) , (74) 



/oo 
E{ai,Si,u;)[l - 6{Si - ai)]B''^''-^\Si,u;)dSi , (75) 
-oo 

from which it becomes apparent that the solution B^^"'\ai,uj), to any order of approximation, 
is expressed as a fraction, the denominator 'D{ai,uj) of which (not depending on the order of 
approximation), can become small for certain values of ai and uj so as to make B^^'^\ai,uj), and 
(possibly) the field in the substratum, large for these values. When this happens, a natural mode 
of the configuration, comprising the blocks, the soft layer and the hard substratum, is excited, 
this taking the form of a resonance with respect to B^^''^\ai,uj), i.e., with respect to a plane wave 
component of the field in the substratum. As B^{ki,uj) is related to A^{ai,ijj) and B^[ai) via 
(|53|) -(|55 |) . the structural resonance manifests itself for the same cti and uo as concerns the field in 
the layer. 

We say that S^(")(cri, w), and the fields in the layer and substratum, can become possibly large 
at resonance because until now we have not taken into account the numerator M^'^\ai,uj), which 
might be small when the denominator is small, or such as to prevent, by other means, the fields in 
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the layer and substratum from becoming large. Moreover, since the field is expressed as a sum of 
plane waves, the fact that B^{ki,uj) may become large for some k^, does not necessarily mean that 
the sum of plane waves (including waves whose horizontal wavenumber ki ^ k^), and therefore the 
field, will be large at a resonance frequency. 

6.2 B^{ki,ij) in the absence of the blocks 

The easiest way to account for the absence of blocks is to take 5 = 0. Then: 



cos{klh) + i-T-^ sin{klh) 



fi^k^ 



D{ki,Ki,uj) = ; VA;i,Ki G 

so that yields 

C{ki,uj) Air 

■ the case of bulk plane wave excitation at ir 
I distribution 6{x — y)F{y) = 6{x — y)F{x), 



(76) 



(77) 



cos{klh) +i^sm{klh) 
cos(A;|/i) - i^sm{klh) 



Vfci G M . (78) 



First consider the case of buU 
of the Dirac delta distribution 6{: 

B"{ki,u;) = — " ^-4- 

47r 



excitation at incidence angle 0*. Using the property 
we obtain 



wherein 



r z 

cos(A;2'*^) — i ^o^i,i sm{k},'^h) 



6{ki - k\) ; VA;i G M , (79) 



4\2 



^ (A:^'') > , 9 (a;^-') > ; i = 0, 1, 2 . (80) 



What this result means is that the amplitude function B^{ki,uj) vanishes for all ki except ki 



k 



1' 



and since (if we assume that there is no dissipation in the layer and substratum), ii^k2'^ / ^i^k^'^ , /cg'*) 
A;2'*, sin(A;2''/i), cos{k\'^h) are real, the denominator in ()79() cannot vanish (neither be small if there 
is a reasonable amount of dissipation in the layer and / or the substratum) , which is another way of 
saying that no structural resonances can exist when the configuration (without blocks) is excited by 
a bulk plane wave. 

In fact, it is easy to ascertain that the existence of structural resonances, in the case of the 
configuration without blocks, is tied up (when there is no dissipation in the layer and substrate) 
with the possibility of /cg becoming pure imaginary, because then the denominator can effectively 
vanish (or become very small for a reasonable amount of dissipation in the layer and/or substratum). 
This possibility is connected with the excitation of a Love mode, characterized by the simultaneous 
existence of a surface wave (associated with pure imaginary ^2) in the substratum, and a standing 
bulk wave (associated with pure real ^2) in the layer. 
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To make this more palpable, we introduce (jTH]) into H39|) so as to find 



u (x, Lo) = S{uj) exp < i klxi 



1 0,i 



+ 



S{u;) 



003(^2'*^) + i ^o^il' sin(A;2'*/i) 



cos(A:2'*^) 



- sin(A;2'*/i) 



exp |i fcixi + k2''ix2 - 2/i)j | ; Vx e 



(81) 



which expresses the fact that the field in f^o is the sum of the incident plane wave and a reflected 
plane wave for all k^. Both of these waves are bulk waves (i.e., both of the cartesian components of 
their wavevectors are real for non-dissipative media) as opposed to the requirement of there being a 
surface wave (i.e., whose X2 component is imaginary for non-dissipative media) in when a Love 
mode is excited. This means that it is impossible to excite a Love mode, in a configuration without 
blocks consisting of a soft layer overlying a hard halfspace, when the incident wave is a plane bulk 
wave. 

However, as shown in |2H I22j. it is possible to excite a Love mode, in a configuration without 
blocks consisting of a soft layer overlying a hard halfspace, when the incident wave is a cylindrical 
wave radiated by a shallow source. 

In the next section, we shall show that it is possible to excite something like a Love mode 
in a configuration consisting of a set of buildings in welded contact with a soft layer overlying a 
hard halfspace even when the incident wave is a plane bulk wave, the same being true when the 
solicitation is a cylindrical wave. 

6.3 Quasi Love modes 

Let us return to the denominator of the expression of which, for convenience, we re-write 

in terms of ki: 



V{ki,u)) = C{ki,uj) - k^D{ki,ki,uj) 



Ui 



cos{klh) — i^-QT§ sm{k2h) — cos{k2h) — i^-j-rf sin{k2h) 



1 

2^ 



m=0 



tMk}£b,)lji^-{ki,cv)I^^+{ki,uj) ■ VfciGM, (82) 



wherein it is easy to show that 



I^^^{k^,u.) = '-e^''^' 



smc I + — ) + (-1) smc f 



rmr\ 



(83) 



with 



sinc(C) := 



sin(C) 

c 



(84) 



Since C{ki,u;) := C{ki,uj) = is the dispersion relation for (i.e., providing the means of deter- 
mining the (ki,ij) couples leading to a possible resonance associated with the excitation of) Love 
modes, we can say that (|H^ is the dispersion relation for quasi Love modes. 
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Remark 1 

Quasi Love modes are different from Love modes wliich, at present, means that the {ki,uj) couples 
for which T>{ki,u) = are not identical to the {ki,ij) couples for which C{ki,uj) = 0. 

Remark 2 

When 6 — > 0, the dispersion relation for quasi Love modes becomes the dispersion relation for Love 
modes. 



Remark 3 

For small k^Q the quasi Love modes are a small perturbation of Love modes, which means here 
that the {ki^oj) couples for which V{ki,uj) = are close to the (fci, couples for which C(A;i, w) = 0. 

Remark 4 

For small fj?/ fj,^, the quasi Love modes are a small perturbation of Love modes, which means here 
that the ki,uj couples for which D^kijUj) = are close to the ki,u; couples for which C{ki,uj) = 0. 

Remark 5 

For small k'^w the quasi Love modes are a small perturbation of Love modes, which means here 
that the ki,io couples for which D^kiju) = are close to the ki,u; couples for which C{ki,io) = 0. 

Remark 6 

The dispersion relation for quasi Love modes is independent of the number of blocks (provided this 
number is greater than 0). 

To substantiate these remarks (when necessary) and obtain a more detailed picture of the fea- 
tures of the quasi Love modes as compared to those of the Love modes, we must analyze more 
closely ()8'2() . For m = we have 

/(^■)^(A:i,..) = e±i'=^^sinc(^) , (85) 

and for kiw ^ zbmvr, we have 

/W±(/ci,^) = ^ ^ ^e±ifci-sm^^ + — . (86) 

[kiWj) — [rmr) V ^ ^ / 

To make a complex issue relatively simple, we assume that kiw is effectively such as to be different 
from rmr for m = 0, 1, 2, .... Then 



'D{ki,uj) = cos^k^h) — i ^Q Q sm{klh) — ^q, q cos^klh) — isin{klh) 

(Ji tXiry LL rXir) 



^2 



3> 



27r Li^ I kk y 2 / ^ ^ ko /t \2 / \2 

r~ \ z \ / z [Kiw) — (rmr) 



. fkiWj rm\:\ ( k\w^ rmx\ ,,2(i), A , 
sin + —J sm (^-^ + —J tMkZ^b,) \ , (87) 
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which rather clearly substantiates the aforementioned remarks. 

Consider the first term in {..}. This term is significant only for small kiWj/2 due to the sine 
function whose modulus decays rapidly as its argument increases. Another feature of this term is 
that it vanishes when k'^bj = Itt ; / = 0, 1, 2, which occurs when the zeroth-order quasi-mode in 

a block encounters a stress-free boundary condition at the base of the block (i.e. u^2'''' — 0)) 

' X2=0 

in which case the latter is disconnected from the underground (since no wave can penetrate into 
the layer) from the point of view of the fundamental block quasi mode. This corresponds to the 
stress-free base block mode. It is then logical that this quasi mode of the block should not perturb 
the dispersion characteristics of the mode of the whole configuration. 

The analysis of the series term in {..} is more difficult. It is clear that a few terms of this 
series should be retained, unless k'^b and/or kiw are very small. The subsequent terms of the series 
become rapidly small (with m) due to the fact that tan(/c2^'*6j) ~ itanh^rmrbj / Wj) i ; m ^ oo 

and T7r^%^T^W = 0{m^^) ; m ^ oo. 

[{kiWj) -{miT) \ 

In any case, it seems legitimate to adopt the following picture of what is going on: the base 
of a given block is a location where diffraction is produced resulting in incident (from either the 

block into the underground or from the underground into the block) bulk waves being transformed 

(i)+ (i)— 

into diffracted bulk and surface waves, as is testified by the presence of the terms Im and Im 
in the series, and by the fact that these diffraction effects disappear as the width of the base 
segment goes to zero. Naturally, the presence of these locally-produced diffracted waves perturbs 
the overall wave structure (with respect to what it was in the absence of the blocks) and therefore 
results in a modification of the characteristics of the modes in the layer and substratum (which 
were Love modes when the blocks were absent). This picture is consistent with the observation 
that the diffracted waves are more difficult to produce when the base segment of a block appears 
as a stress- free surface due to the fact that either k'^bj = Itt ; / = 0, 1, 2, ... or ^^///^ << 1. 

Beyond this, it is necessary to carry out a numerical study in order to see how the different 
parameters involved in the problem, notably those of the block, modify the dispersion characteristics 
of the modes of the configuration with respect to what these modes were (i.e.. Love modes) in the 
absence of the blocks. The numerical study should also seek to evaluate the modification of the 
response (notably on the ground) of the configuration due to the presence of blocks. 

6.4 The emergence of the natural modes of the configuration from the hnear 
system of equations for Bm''^: quasi displacement-free base block modes 

Eas.((7T |) - ((7^ signify that ^i?^^'^^ becomes large when C„ — Dnn becomes small, and that 

this occurs at all orders p of approximation. The fact that ( Bn j becomes inordinately large 
is associated with the excitation of a natural mode of the configuration. The equation Cn{uj) — 
Dnn {(^) = is the dispersion relation of the n-th natural mode of the configuration. 
Let us examine this relation in more detail. 



-<k'^S) - '-^^i^ j j^^Hk,,.)i^^-{k,,.) 



cos{k^h) - i^sm{k^h) 



cos{klh) - i^sin(A:^/i) 



11,1 



dki 

To" 

«,2 
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which shows that the natural modes of the configuration result from the interaction of the fields in 
two substructures: the superstructure (i.e., block(s) above the ground), associated with the term 

.FS(6,/i^p^u;)=.FS=cot(4')60 , 



and the substructure (i.e., soft layer plus hard half space below the ground), associated with the 
term 



-r-Wfu 1 1 2 2 \ T-i 



2n 



■ 2; 2(0 
27r/i0 



/«+(fci,u;)/«-(A:i,a;) 



ms{klh) - i^sin(A;|/i) 



cos{klh) - \^sm{klh) 



dki 
To" 



(90) 



Each of these two substructures possesses its own natural modes, i.e., arising from Tin{b, /x^, p^, w) = 
for the superstructure, and J^2n{h, , , , , fi^ , ,w,(jj) = for the substructure, but the 
natural modes of the complete structure ( superstructure plus substructure ) are neither the modes of 
the superstructure nor those of the substructure since they are defined by 



.fS (6, /i^ p^ a.) - .F^;^ (/i, p\ /x^ p^ T., u;) = , 



r(0 



„o 



(91) 



which emphasizes the fact that the natural modes of the complete structure result from the inter- 
action of the natural modes of the superstructure with those of the substructure. Note that the 
solutions of (|9T|) can be quite different from those of either /both J^in = and J^2n = 0. 

In order to establish the natural frequencies of the complete structure, we first analyze the nat- 
ural frequencies of each substructure and assume that all media in the structure are non-dissipative 
(i.e., elastic). 

The solutions of the dispersion relation for the superstructure are: 



cot 



4^ k^n^l = (2m- + 1) TT 44> UJ = UJn^ = C" 



bi J \wi 
; m,n = 0,1,2,... , (92) 



which are the natural frequencies of vibration of a block with a displacement-free base (i.e., at these 
natural frequencies, ?x^'''^|^^_q vanishes on the base segment of the block). This relation shows that 
there corresponds to the n-th block mode an infinite set of sub-modes, identified by the index m. 
The same can be said about the natural modes of the entire configuration. 

Now consider the dispersion relation for the substructure of the entire structure. Due to the 
fact that the integrand is an even function of ki, it becomes 



■ 2 7 2(0 

(0 _ k2„'wi 



2n 



• 2; 2(0 



/«+(A;i,a;)/W-(fc,,^) 



/«+(A;i,u;)/«-(A:i. 



cos 



sm 



(klh) 



cos(A;2^) 



sin{k2h) 



dki 
To" 

"-2 



p^k2 cos{k2h) — ip^k2 sm{k2h) 
p^k^ cos{k2h) — ip^k2 sm{k2h) 



dki 

^2 



(93) 
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A few preliminary remarks are in order: 



Remark 1 

Since the term in the dispersion relation is absent in the absence of the infrastructure, we can 
say that the natural modes of the complete configuration are quasi displacement-free base block 
modes. Quasi displacement-free base block modes are different from displacement-free base block 
modes which, at present, means that the {n,u)) couples for which Tin{b, fx^ , ,uj) = are not 
identical to the (n,a;) couples for which J^2n{h, jjP , , fJ-^ , , iJp , ,10,0;) = 0. 

Remark 2 

For small p? / p}, the quasi displacement- free base block modes are a small perturbation of displace- 
ment- free base block modes, which means here that the (n, w) couples for which !Fin{b, p'^, p'^,u>) = 
are close to the (n, co) couples for which T2n{h, p^, p^, p^ ,p^ , p'^ , p"^ ^w, w) = 0. This is a relatively- 
logical result in that when p^ / p^ is small, the waves coming from the infrastructure have more 
trouble penetrating into the blocks and modifying therein the modal structure. 



Remark 3 

The dispersion relation for quasi displacement- free base block modes is independent of the number 
of blocks (provided this number is greater than 0). This is a somewhat surprising result related 
to the choice of the iteration method for solving the linear system of equations for B n , Since a 
more accurate choice of method (one of which is described in sect. 16.5(1 can be shown to lead to 
a somewhat different (although much more complicated) dispersion relation which depends on the 
number of blocks in the configuration. 



We now analyze in more detail J^2ni and, in particular, J^20- Recalling (jH^ . we get: 



P 

■J 



(I) _ ik'^wi p"^ 



20 



TT p' 



smc 



cos(/c|/i) - ijj^sm{klh) 
cos{klh) -i^sm{klh) 



dki 



Proceeding as |211 122j . we decompose the integral into three parts (under the assumption 
> 0) so as to obtain: 



-rW/j, 1 1 2 2 N T-^(l-) , -r-2(«j I ^■ 

•^2o(^'/" '/^ 'P ^^'^) =-^20 +-^20 +-^:: 



1(0 



-.2(0 



-3(/) 
20 



wherein 



1(1) _ iPw p^ 



20 



vr p 



T' 



.2(0 _ k'^w p^ 



20 



smc 



smc 



T' 



3(0 _ k^w p^ 



20 



vr pi 



smc 



cos(Ki/i)-i^sin(i(-i/i) 



co^{K\K) - ij^ w^(K\h) 

"cos(Ki/i) + ^sin(A:i/i) 
co^{k\K) - |^sin(A;^/i) 



dk\ 



cosh(K2^) + ^^j-j^ sinh(K2/i) 
cosh(K2/i) + ^"0^ sinh(K2/i) 



dki 



(94) 
k^ > 
(95) 

(96) 
(97) 
(98) 
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with 

kiioj) = Kiico) := ^{ki{oj)f - {kif ; ki < k{{uj) ; w > , (99) 



ki{u) 



{k,f-{k^{oo)y 



ki > k^{uj) 



uj> 0. 



(100) 



As shown in [211 \'22\ . usually dominates the other two terms, and this is due to the fact that 
the denominator in the integrand of can vanish for ki over a large portion of the interval of 
integration for frequencies at which Love modes are excited in the infrastructure, this occurring 
near the Haskell frequencies 

,,LOVE _ ,,HASK _ + 1 



'I 



2h 



m G N , 



which corresponds to 



■.0 LOVE 



■0 MASK 



Ittu^"^^^ _ 2m + 1 TTC^ 
^0 ~ 2 



m G N . 



(101) 



(102) 



k 



1 LOVE 



f^l HASK 



HASK 



2m + 1 vr 

; mGN 



(103) 



The sinc^ function in the integrand of J^2o is significantly large only in the interval [0,27r/w], so 
that a minimal requirement for capturing most of the contribution of the Love modes in ^lo 
their frequencies of resonance is that 

2tt 



w 



(104) 



Of course, there exist other terms in the dispersion relation of the n = mode, i.e., J^20j -^Io' 



and the contributions of the higher-than-fundamental Love modes to .T-fi 



As concerns J^^n and 



2Q. JT-O V^Uilv^CillD ^ 20 

we notice that the former is complex and the latter is pure imaginary, whereas J'Iq, is real, so 
that we would expect to have J-2Q and J^|q to contribute less to the dispersion relation than J^2o- 
The contribution to the n — th order natural mode of higher-than-fundamental Love modes to T^q 
is empirically found to be always less than that of the zeroth order Love mode. 

The analysis of the n > natural modes of the complete structure proceeds in the same manner 
as previously, and will not be given here. 



6.5 Another look at quasi displacement-free base block modes 

The system of linear equations (|66j) can be written as 



wherein 



2er 



cosiklilhi) 



1 m=0 



4''+(fci,c^)e'^^('^'-""/^) 
cos(/c^/i) - ij^isinik^h) 

; V/ G B ; Vn G Z , (106) 



dki 
To" 

"'2 



(105) 
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and 



21.2 



2m 



2^ cos{kl!;^hj) 

/oo 
/«+(fei,a;)jC^)-(A;i,u;) 
-oo 



cos(klh) — i^Vrr sinf/cn/i) ..lui 



cos(A;2^) 



; V/, j G B ; Vn, m G Z . (107) 



Remark 1 



If = dj, then Q^il = for m + n = 2p + 1, Vp G N. 

Remark 2 

If = 0, then = forn = 2p + 1, Vp G N. 



Remark 3 

2(1) 

According to remarks 1 and 2, -82^+1 = 0, Vp G N, when the source is located on the vertical line 
passing through the center of the base of the lib. bloc. This block then only admits displacement 
via the even quasi-modes. 



Let M := {1, 2, 3, A^} and 

b:= (b^W b2(2).... h^Wy , 
where T is the transpose operator, and 



Similarly, let 



b2(0 (^^2(0 ^2(0 ^2(0 ____ Y 

p:= (p« p«.... pW)^ , 



p(0 (^pW p(v p 



(0 d(0 

2 



(108) 

(109) 

(110) 
(111) 



Q:= 



Q(liV)x 



Q(21) Q(22) Q(2JV) 



\^Q(iVl) Q(iV2) ____ Q{NN)J 



(112) 
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Q 




qW 


Q 


^20 


QIC 


Q 



Q 



02 
(Ij) 
12 

22 



V 



(113) 



so that pU5|l can be written as 

b = p + Qb, (114) 

This is a matrix equation enabUng the determination of the set b of unknown vectors {b^^^^ ; / € 
B}; it can be written as 

(I-Q)b = p, (115) 

wherein I is the identity matrix (i.e., the matrix having the same dimensions as Q with one's on 
the diagonal and zeros elsewhere). 

The natural modes of the city are obtained by turning off the excitation, embodied in the vector 
p. Thus: 

(I-Q)b = 0, (116) 
which possesses a non-trivial solution only if 

det(I - Q) = , (117) 

wherein det(M) signifies the determinant of the matrix M. 

To get a grip on this relation, consider the case in which there is only one block in the city. 
Then (|117|) becomes 



det 



V20 



V 



Vol 
— Q21 



^02 
— Ql2^^ 
1 ~ Q22 



, 



(118) 



A procedure, called the partition method, for solving this equation (as well as the equation (|115j) for 
the response) , particularly appropriate if the off-diagonal elements of the matrix are small compared 
to the diagonal elements, is first to consider the matrix to have one row and one column. 



then to consider the matrix to have two rows and two columns, 

'1 



det 



VlO 



1 



-Qoi ^ 







(119) 



(120) 
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or 



(11)^(11) 



Qni Q 



10 



then to consider the matrix to have three rows and three columns, 



:i2ii 



det 



V -Q 



00 

Vio 

(11) 

20 



Vol 
-Q21 



"Q12 
~ Q22 



:i22) 



or 



Voo 



^(11)^(11)^(11) 
Vio V21 V02 



Q22 



^(11)^(11)^(11) 

V20 V12 Vol 



Qn I Q20 ^Qo2 ^ 



Q22 I Qio ^Qoi ^ 



0, (123) 



etc. 

If the off-diagonal elements of the matrix are considered to be negligible compared to the 
diagonal elements, the dispersion relation 1)118^ is simply that the product of the diagonal elements 
of the matrix should be nil, which means that any diagonal element of the matrix should be nil. 

Consider the case in which it is 1 — Qqq^ that vanishes. This is equivalent to 



cot{k'^b) 



• ; 2 2 

IK Wi H 

27r 



cosik^h) - i^ sin(A;^/i) 



cosik^h) - ij^sin(fe|/i) 



dki 
To" 



(124) 



which is the same as the dispersion relation of the zeroth-order quasi displacement-free base block 
mode obtained in the previous section. Thus, there is no apparent gain in adopting the analysis of 
the present section over that of the previous section in an attempt to resolve the difficulty mentioned 
in Remark 3 in the previous section, namely the fact that the dispersion relations do not depend 
on the number of blocks in the city. 

To resolve this problem, we must take into account the off-diagonal elements of the matrix be- 
cause these elements contain the information on the number of blocks. Unfortunately, the inclusion 
of these off-diagonal elements makes the dispersion relation increasingly complicated and difficult to 
analyze as the order of approximation of the partition method (which consists in solving increases 
(it is already quite complicated at first order). The only way to solve these dispersion relations 
(which consist in equating determinants of increasing rank to zero) is by numerical computation. 

When more than one block (for example 2 blocks) are present, the dispersion relations becomes 
(if the off-diagonal elements of the matrix for each block considered independently can be neglected): 



det 



1 



Voo 



Voo 



Vqo 



n(22) 

Vqo 



1 



Voo 



1 



^(12)^(21) 

Voo Voo 







(125) 



which is quite different from the zeroth-order quasi displacement-free base block dispersion relation 

(12) (21) 

because the coupling term Qqq Qqq does not vanish and cannot be neglected. This term corre- 
sponds to the so-called structure- soil- structure interaction and accounts for the distance separating 
the two buildings. Its form is close to that of the term representative of the action of the geophysical 
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structure. Nevertheless, due to its complexity, it is difficult to carry out an analytical study of this 
relation. 

The partition method emphasizes the role of the global superstructure, while the iteration 
method, described in sect. 16.41 emphasizes the role of only one component (i.e. one block) of 
the superstructure. The partition method also accounts for all the possible interactions between 
blocks (the term cot (/C20&O) ^-^d the geophysical structure through the terms Qnm, and the inter- 
action between blocks through the terms Qnm and Qiim- 

Ultimately, the choice of method reduces to determining which one gives the best results, i.e., 
results that are closest to reality. This can be determined only by full-blown numerical studies, the 
results of which will have to be compared to those of the FE method (employed as a reference). In 
fact, we find that the partition method gives the best results, and is therefore employed in all the 
subsequent numerical computations. 

7 Expression of the fields w^'^^)(x, cj), w^(x, cj) and w^(x,a;) for line 
source excitation 

Once the quasi-modal coefficients Bm , '^rn G N, V/ G B are obtained from the matrix equation 
the field in the block domain O^^'^ is computed via (|44|) . This field vanishes on the ground 
at the frequency of occurrence of the displacement- free base mode of the block. 

Remark 1 

If higher-than-the- zeroth-order quasi modal coefficients can be neglected, the field in the block 
/ G IB takes the form 



u 



2(0 ^ o2(0 



b'^^'>{uo) COS {k^{x2 + hi)) (126) 

which indicates that the displacement field is independent of xi and takes the form of a standing 
wave in the block. 

Let us next consider the field in the layer. Combining (|T7|) . (|5T|) . (|53|) and (|55|) . leads, via (|ifl|) . 

to: 

n"^(x, uj) = 

-.Si.) r - (*■!-) .'<'--»«»'-^-") i E E Sin (4!;'.) X 



2» cos(tifc)-i^dn(tlft) ^'»if„f° 



COS {k\ix2 — h)) + i^^^r^ Sin ikl(x2 — h)) . jl^ 

^ , [ \ — ^/«-(^„.)e-^(-(^-^))^ , (127) 

cos {k\h) -\^^s\n{klh) 
which can be cast into the form 

u\^,io) = ulix,u;)+Y,UM\^,^) , (128) 

with 

1, , i5(a;) Z''- cos (A;ix2) e^(''^("^-"i)+'2("2-'^)) , , 



2^ ■'-'^ cos (Ajg/i) - i^^ sin (A;^/i) ^2 
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and 

oo 2 
n=0 ^ 



1(0 r \ 



cos(A;i/i) -i^sin(A;i/i) ^^2 



This expression indicates that the field in the layer is composed of: i) the field obtained in the 
absence of blocks and induced by the incident cylindrical wave radiated, ii) the fields induced by 
the presence of each block. The displacement field n^(x,a;) — uj(x, w) appears as a sum of block 
fields which are strongly linked together, since each coefficient Bn is calculated by taking into 
account the presence of the other blocks via pi5|) . 

Let us finally consider the field in the substratum. Combining (|17|) . (|5T|) . (|53|) and (j55|l . leads, 
via dnni), to: 

nO(x,u;) = nXx,a;) + T ^^^^^ + ^^^^^ ^i(,,(.,_,..)+fe0(,^+,._2,)) dfci ^ 

4^ J~oo cos rfci/v.^ - «in rfclM ^2 



A;i/i)-i^sin(A;i/i) ^^2 



00 



«6 

which can be cast into the form 



„=o ^ ^ ^ ioo COS (A:^/i) - 1^ sm (k^h) '^2 



^i°(x,w) = {u^(x,w) +nf (x,u;)} +^nj'^(x,cj) (132) 

with 



i^(^) cos (fc^/.) + ig| sin (fc2^fe) ^.^^^ _.,,^,0(.,^.,_2.)) rffc 
4^ ^-oocos(A:i/i) -i^sin(fci/i) ^2 



and 



nr(x,..) = ^E;^^^«^?..,sm(4'^6. 



n=0 



cos(fc^/i) -ij^sin(/ci/i) ^2 



(134) 



This expression indicates that the field in the substratum is composed of: i) the field obtained 

in absence of blocks, including the incident plus diffracted fields, the latter being induced by the 

incident cylindrical wave, ii) the fields induced by the presence of each block; the latter takes the 

2(7) 

form of a sum of block fields which are strongly linked together, since each coefficient Bn is 
calculated by taking into account the presence of the other blocks via (|115() . 
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7.1 Interpretation of the fields u^''^ and u^''^ 

If the leading term of the quasi modal representation in each block / G B is dominant (i.e., the 
higher-order terms can be neglected), the fields u^^\'k,uj) in (|13()|) and u^^\'x.,io) in (|134j) reduce 
to: 



u 



1(7). ^ i 2(i)r2 . 2, N /■°°sinc(A;i^)ei('=i(^i-^:')+'^2(^2-/i)))^^ 

2vr/.o ^ ^ ^^ioo cos(fci/.)-ig|sin(A:i/i) ^2° 

(x, = ^Bq- 'k Wjsm[k bj) / — -g- . (136) 

2^/^ Joo cos (klh) - i sin (klh) ^2 

With the help of the material in appendix (^, each v!^''\'x.,uj) can be interpreted as the field 
radiated by ribbon source of width Wj, located at the base of each block j G B, and whose amplitude 

is of the form ^B'^^''^k'^)wj sin {k^hj) . These sources are induced sources, i.e., they do not introduce 

energy into the system, but each of them induces a modification of the repartition of the energy 
over the excitation frequency bandwidth. These sources are located at the top of the layer and 
should excite quasi-Love waves, as shown in [221221120] (for applied sources). 

If higher-order quasi modes are relevant, the even-order modes corresponds to ribbon sources 
of width wi, located at d/, / G B, while the odd-order modes correspond to line sources located at 
the edges of the blocks. The amplitudes of both of these types of induced sources depend on the 
order of the quasi modes and on the characteristics of the corresponding block. 



8 Numerical results for one block in a Mexico City-like site 

All intervals over the which integration is performed in the calculation of the quasi-modal 
coefficients, through the linear system (|115|) . i.e. in the calculation of Pn"* and of Qnm^ {n, m) G N^, 
(j, /) G B^, and in the calculation of the displacement fields = u\ + ^^w^'^ and vP = + 

P + OO 

u^'^ + 5^ % ' reduced to / . These intervals are then separated into Ii = [0, [k^)] 

zeB 

(interference of propagative waves), X2 = {k^) ,^ [k^]] (excitation of quasi-Love waves) and 
^3 = {k^) 1 +00] (evanescent waves in the layer) in order to point out the different type of waves 
associated with the different possible phenomena in the geophysical structure. The numerical 
evaluation of these integrals is carried out by the the procedure described in j22] . 

The results in this section apply to a single block whose base segment center is located at 
(Om,Om). 

We compute the response inside the block and on the ground near the block. 

The latter is supposed to be situated in a Mexico City- like site wherein: = 2000 kg/m^, 
€"=600 m/s, = 00, = 1300 kg/m^, c^=60 m/s, = 30, with the soft layer thickness being 
h = 50m. In addition, the material constants of the block are: p = 325 kg/m^, c2 = 100m/s, 

g2 = 100. 
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Figure 6: Indications concerning the solution of the dispersion relation T = T\— J-2 = for a single 
40m X 40m block in a Mexico City- like site. In the left panel the solid curve describes (^i) versus 
frequency {y = u;/27r in Hz) and the dashed curve describes 5? (^2) versus frequency. The right 
panel again contains two curves: the solid one describes 9 i^i) versus frequency and the dashed 
curve describes 9 versus frequency (Hz). 



The source is placed consecutively at (Om, 3000m) or (-65m, 3000m), which are deep locations 
for which Love modes can hardly be excited in the absence of the block, and at (-3000m, 100m), 
a shallow location at the which Loves modes can easily be excited in the absence of the block. 



The eigenfrequencies of the block displacement-free base block are z^om 

2m + 1 



_ c2(2m + l) 



2b 



and the 



Haskell frequencies are 



HASK 
m 



wherein m = 0,l,2,.... Thus, the Haskell frequencies 



2 2/1' 

are 0.3, 0.9, 1.5 ...Hz. The fundamental displacement-free base block eigenfrequency (whose value 
is supposed to be close to the one of the corresponding quasi-mode), depends on the choice of the 

25 

— ~ — . This expression 
2o 



mechanical parameters of the medium filling 0.2, and occurs at y^^f^ 



agrees with the empirical one: 



FB 

00 



30 



employed in |^. 



If the zeroth-order quasi-mode is dominant, the dispersion relation (|124)) takes the form: 



= cot(A;^6) 



2m _ 

27r 



cos{klh) - ij^sin(A;|/i) 
cos{klh) - ij^sin(A;|/i) 



2 / , 'w\ dki 
smc ih^j -j^ = ^1- ^2 



(137) 



8.1 Results relative to one 40m x 40m block 

The block is 40m high and 40m wide. The displacement-free base block eigenfrequencies are 
then 0.625, 1.875 Hz.... 

Fig. ini gives indications on the frequencies of occurrence of the modes of the entire (superstruc- 
ture + substructure) built site. An eigenfrequency is a frequency for which = Ti — T2 = 0, i.e., 
for which ?R.{J^i) = ?R.{J^2) at the least. One notes that this occurs for ~ 0.3, 0.63, 1.35, 1.92 Hz 
in the frequency range of the figure. An indication of the attenuation associated with a particu- 
lar mode (at a frequency i^*) is the quantity |9(J^i(z/*)) — Q {T'2{'^*)) \ - The quasi-Love mode at 
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27rI|cos(k^b)x(1-QJ|| 

2rr||cos(k^b}|| 

27r| cos(k^b)xQ|^i^|[, integration over k^E 1^ 

27r| cos(k^b}xQ^^||, integration over kjE 1^ 

_ 2j[||cos(k^b}xQ||j|||, integration over k^e 1^ 
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Figure 7: 27r|| cos(/c^6)(l — (5oo)|| (left panel) and the driving agent 27r|| cos(A;^6)Po/'S'(w)|| (right 
panel) . 

u ~ 0.3Hz and ~ 0.9Hz, and the quasi displacement-free base mode block at i/Q^^^ 0.63Hz 
and ~ 1.9Hz are modes with relatively-small attenuation. On the contrary, the quasi stress-free 
base block (close to stress-free base block modes, corresponding to M^^^g = and satisfying 

tan (/c^6) = 0) at I'^i^^ ~ 1.3bHz is a mode with very large attenuation. This third type of mode 
(quasi stress-free base block modes) should therefore have a small influence on the response at the 
site, as was mentioned in section 

To give an indication on how the block is excited, and how it re-radiates the field, we depict 
in figure [71 the absolute values of the zeroth-order terms, i.e. (1 — Qoo) and Pq; involved in the 
resolution of the linear system (|115j) . The integration interval is divided into three subintervals: 
Xfc, k = 1,2,3. The source term 27r|| cos(/c^6)Po('^)/5'(i^)|| is similar to the transfer function in the 
absence of the block in |22j . The term 1^^ {ki,io) = sine (^i-f ) exerts a small influence on Pq{uj). 
The solicitation of the block takes a form close to waves traveling in the layer in the absence of the 
block. The deep source being located close to a line going straight down from the block, the main 
component of Po{uj) comes from the integration over Xi (i.e. interference of propagative waves 
in the substructure), while the main component of Qoo{^^), which is independent of the source 
location, comes from the integration over T2 (i.e. Love mode excitation). 

8.1.1 Displacement field on the top and bottom segments of the block for deep line 
source solicitation 

We now examine the displacement field on the horizontal boundaries of the block. 

Figure |H1 depicts the spectrum and time history of the total displacement at the center of the 
summit segment of the block as computed by the mode-matching method (with account taken of 
one quasi-mode) and the finite-element method for a deep source. No noticeable differences are 
found between the results of the two methods of computation. The neglect of the quasi-modes of 
order larger than is valid for this block width. The block acts as a ribbon source of width w 
located at the base segment. 

Figure depicts the spectrum and the time history of the total displacement on the ground in 
the absence of the block as well as the total displacement at the centers of the top and base segments 
of the block for a deep source. In the time domain, a small increase of the duration, and a fairly- 
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■ V (Hz) ■ ■ ■ • ,(s) 

Figure 8: 2tt times the spectrum (left panel) and time history (right panel) of the total displacement 
response to the cylindrical wave (radiated by a deep source located at = (—65m, 3000m)) at 
the center of the summit segment of a single 40m, x 40m block. The dashed curves correspond 
to the semi-analytical (mode-matching, one mode) result and the solid curves to the numerical 
(finite-element) result. 




Figure 9: Comparison of 2it times the spectrum (top panels) and time history (bottom panels) of 

the total displacement on the ground in absence of block (solid curves) with the total displacement 
(dotted curves) at the center of the top segment (left panel) and at the center of the base segment 
(right panel) of a single 40m x 40m block, for a deep source located at = (—65m, 3000m). 
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large increase of the peak amplitude, can be noticed, particularly on the top segment of the block, 
where the quasi displacement-free block mode makes itself felt. In the frequency domain, there is an 
increase of the maximum amplitude and the sharpness of the first resonance peak of the substructure 
at around 0.3Hz, as noticed previously, to some extent, in [321 llUj . and is a characteristic feature 
of the so-called soil- structure interaction. A similar phenomenon was mentioned in 16.31 and is due 
to the ability of shallow sources to excite Love modes in the configuration without blocks j21j . 
Here, a quasi-Love mode can be excited due to the presence of the block, which acts like a shallow 
(actually located on the ground) source. A similar fact was already noticed in [221, where the spectra 
of the displacement on the ground for a deep source (provoking only interference phenomena) was 
compared to that for a shallow source (provoking Love mode excitation). 

The so-called soil- structure interaction consists (in the configuration with blocks) in a modifica- 
tion of the phenomenon from a state where interference (when a soft layer is present) or reflection 
(when the soft layer is absent) phenomena are dominant to a state where a quasi-Love mode is 
excited. This mode is excited because of the presence of the block in a configuration. In the absence 
of the latter, no mode can be excited. 

Another point of view is that the block takes the form of an induced source at the top of the 
layer (see section 17. 1() , and it is now known that Love modes are excited when such a configuration 
is solicited by a line source near the boundaries of the layer |2nj . 

A particular feature of the spectrum of the displacement when the block is present, (top-right 

subfigure of fig. IH} is that it vanishes at the base segment for Vqq^^ ~ — r , which is the fundamental 

displacement-free base block eigenfrequency. This is made evident in the field representation in 
the block and can be understood to be a geometrical effect which, in the case of a single block, 
corresponds to the excitation of a quasi-mode. 

8.2 Results relative to one 50m x 30m block 

The block is 50m high and 30m wide. The displacement-free base block eigenfrequencies 

are then 0.5, 1.5Hz, Figure 11171 gives an indication of the modes of the configuration. One 

notes that the eigenfrequencies (frequencies at which 5R(.Fi) = 'Si{J^2), at the least) are v « 
0.3, 0.5, 1.0, 1.53, 1.75, 1.93. The attenuations of the quasi-Love mode at ~ 0.3 and 0.9 Hz, 
and of the quasi displacement-free base block at ~ 0.5 and l.bHz are relatively-small. As 
previously, the attenuation of the quasi stress-free base block (close to the zeros of tan [k'^bj) at 

^01^^ ~ \Hz is relatively-large. 

In figs. Illl andlT^ l we plot the absolute values of the zeroth-order terms involved in the resolution 
of the linear system H115|) . The source terms 27r|| cos(A;^6)Po('^)/'S'(w)|| are once again close to the 
transfer functions calculated without the block for both locations of the source = (Om, 3000m) 
and = (-3000m, 100m). 

The term 1^^ (fci, w) = sine (^i^) again has a small influence on Pq{uj). For a deep source, the 
main component of P{){uj) comes from the integration over Zi (i.e. the solicitation of the block takes 
the form of interfering propagative waves traveling in the layer associated with bulk waves in the 
substratum) while for a shallow source, the main component of ^0(1^) comes from the integration 
over I2 (i.e. the solicitation of the block takes the form of propagative waves traveling in the layer 
associated with evanescent waves in the substratum, i.e.. Love modes). The integral over I2 again 
dominates Qoo('^)- 
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Figure 10: Indications concerning the solution of the dispersion relation ^ = — ^2 = for a 

single 50m x 30m block in a Mexico City- like site. In the left panel, the solid curve describes 3? (^1) 
versus frequency {i' in Hz) and the dashed curve describes ^ {^2) versus frequency. The right panel 
again contains two curves: the solid one describes 9 (jTi) versus frequency and the dashed curve 
describes versus frequency . 
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Figure 12: Driving agent 27r|| cos(A;^6)Po/'S'(i^) 
panel) and for a shallow source located at 



for a deep source located at 
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Figure 13: 2tt times the spectra (left panels) and time histories (right panels) of the total displace- 
ment response to a cylindrical wave (radiated by a deep source located at = (Om, 3000m)) at the 
center of the summit segment (top panels) and at the center of the base segment (bottom panels) of 
a single 50m x 30m block. The dashed curves correspond to the semi-analytical (mode-matching, 
one mode) result, and the solid curves to the numerical (finite-element) result. 



8.2.1 Displacement field on the top and bottom segments of the block for deep line 
source solicitation 

We now examine the displacement field on the horizontal boundaries of the block. 

We first consider that the seismic disturbance is delivered to the site by a deep line source 
located at x^ = (0,3000m). This means that in absence of the block, the displacement field is 
mainly composed of propagative waves in the substratum and interfering propagative waves in the 
layer. 

Figure [T^ depicts the spectra and time histories of the total displacement at the center of the top 
and bottom segments of the block, as computed by the mode-matching method (with account taken 
of one quasi-mode) and the finite-element method, for a deep source. No noticeable differences are 
found between the results of the two methods of computation. The neglect of the quasi-modes of 
order larger than is valid for this block width. The block acts as a ribbon source of width w 
located at the base segment. 

Figure depicts the spectra and the time histories of the total displacement on the ground in 
the absence of the block as well as the total displacement at the centers of the top and base segments 
of the block for a deep source. In the time domain, a small increase of the duration (although 
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Figure 14: Comparison of 2tt times the spectrum (on the top) and time history (on the bottom) 
of the total displacement on the ground in absence of the block (solid curves) with the total 
displacement (dotted curves) at the center of the top segment (left panel) and at the center of 
the base segment (right panel) of a single 50m x 30m block, for a deep source located at = 
(Om, 3000m). 
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larger than in the 40m x 40m block case), and a fairly-large increase of the peak amplitude, can 
be noticed, particularly on the top segment of the block, where the quasi displacement-free block 
mode makes itself felt. 

In the frequency domain, there is an increase of the maximum amplitude and of the sharpness 
of the first resonance peak of the substructure, which is a characteristic feature of the soil- structure 
interaction. Again, the displacement field vanishes at the base segment for a frequency close to 

DF 

^00 ■ 

8.2.2 Displacement field on the top and bottom segments of the block for shallow 
line source solicitation 

We now consider what happens when the seismic disturbance is delivered to the site by a shallow 
line source located at = (—3000m, 100m). This means that, in the absence of the block, the 
displacement field in the substructure is that of Love modes at the resonance frequencies of these 
modes. 

Figure El compares the displacement on the ground in the absence of the block to that in 
the block (on the top and bottom segments thereof) for a shallow source. In the time domain, 
the durations are substantially the same when the block is present or absent, but the peak and 
cumulative amplitudes are larger in the presence of the block. In the frequency domain, both the 
sharpness and amplitude of the first peak, corresponding to the first Love mode in absence of the 
blocks, increase. This means that the soil-structure interaction obtained for a deep source is also 
found for a shallow source. The fact that the position of this peak is hardly shifted means that 
the structure of the quasi-Love mode is very nearly that of the Love mode existing in the absence 
of the block. However, the presence of the block enables this mode to be excited more efficiently 
than in its absence, i.e., when the site is solicited solely by the shallow source and not by waves 
re-emitted by the block. 

8.2.3 Displacement on the ground on one side of the block for deep line source 
solicitation 

To obtain a more complete picture of the modification of the phenomena due to the presence 
of a block, we now focus our attention on the ground motion outside of the block. 

The deep seismic line source is located at = (Om, 3000m). Figures ITHl represents 2tt times the 
transfer function (x, to) / S (lj) at various locations on the ground outside of the block, computed 
by the mode-matching method by taking into account only the zeroth-order quasi-mode. From the 
fact that the deep source is close to the vertical line passing through these ground locations, in the 
absence of the block, the first peak of the transfer function is principally due to the integral over 
Ii [21\ \2'2\. The presence of the block induces: i) a double dependency of the first peak on the 
integrals over Ii and over I2 which leads to splitting of the first peak at some locations (as seen in 
the results at Xg = (150m, Om)), ii) surface waves in the layer (represented by the contribution of 
the integration over T3) which normally are not excited in the case of a deep source, |211 \22\ . The 
presence of the block allows the excitation of a quasi-Love mode as is testified by the importance 
of the integration over T2. 

Fig. El shows once again that there is a good agreement between the results of the mode- 
matching method (with only the fundamental quasi-mode taken into account) and the finite element 
method. 
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Figure 15: Comparison of 2tt times the spectrum (on the top), and time history (on the bottom) 
of the total displacement on the ground in the absence of the block (solid curves) with the total 
displacement, (dotted curves) at the center of the top segment (left panel) and of the bottom 
segment (right panel) of a single 50m x 30m block, for a shallow line source located at = 
(-3000m, 100m). 
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Figure 16: 2tt times the transfer function (solid curves) of the total displacement on the ground 
at X = (20m, Om) (top-left), x = (30m, Om) (top-right), x = (150m, Om) (bottom- left), and x = 
(300m, Om) (bottom-right). The deep line source is located at = (Om, 3000m). The dotted 
curve represents the contribution along the interval Ti , the dashed curve the contribution along the 
interval I2, and the dot-dashed curve the contribution along the interval I3. 
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Fig. El shows that when the block is present, the displacement field at a given location on the 
ground is somewhat different from what it is in absence of the block, even at a large (i.e. 300m) 
distance from the block, this being particularly evident in the response spectra. The duration and 
cumulative amplitude of the displacement fields have a tendency of increasing as one approaches 
the block, but not in monotonic manner, as manifested by what happens at 150m. 

This may be one of the causes of the spatial variability of the destructions noticed during 
the Michoacan earthquake. As shown by the form of the matrix in section 16. 5[ one can expect 
one building to have an effect on its neigbors via the waves traveling in the substructure and in 
particular in the layer and along the ground. Thus, the effect of one block on the other can vary in 
a way that depends on the distance of their separation. Nethertheless, the authors of 19 did not 
find substantial differences in response for varying building separations in the case of an idealized 
city composed of ten blocks. It may be that in order for this effect to be noticeable requires either 
a large distance of separation between blocks or a single block. 
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Figure 17: 27r times the spectrum (left panel) and time history (right panel) of the total displace- 
ment at various locations on the ground: (from the top to bottom) at 20m, 30m, 150, and 300m 
from the center of a single 50m x 30m block. The dashed curves correspond to the semi-analytical 
(mode-matching) result (with only the fundamental quasi- mode taken into account), and the solid 
curves to the finite element result. 
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Figure 18: Comparison of 2n times the spectrum (left panel) and of the time history (right panel) of 
the total displacement on the ground in the absence of the block (sohd curves) and in the presence 
of the block (dashed curves) at various locations on the ground: (from the top to the bottom) at 
20m, 30m, 150, and 300m from the center of a single 50m x 30m block. The deep line source is 
located at = (0,3000m). 
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Figure 19: 27r times the transfer function, solid curve, of the total displacement on the ground at 
Xg = (150m, Om) (left panel), and Xg = (300m, Om) (right panel). The shallow source is located 
at X® = (—3000m, 100m). The dotted curve represents the contribution along the interval Ti, 
the dashed curve represents the contribution along the interval T2, and the dotted-dashed curve 
represents the contribution along the interval T3. 



8.2.4 Displacement on the ground on one side of the block for shallow line source 
solicitation 

When the seismic disturbance is delivered to the site by the wave radiated from a shallow line 
source located at x® = (—3000m, 100m), it should be recalled that, in the absence of the block, the 
first peak (associated with the first Love mode) is mainly composed of the component coming from 
the integration over X2, whereas we observe in fig. ^Jthat when the block is present, this peak is 
composed mainly of two components, one from integration over 2i and the other from integration 
over T2. However, this double dependence of the first peak finds no apparent translation in the 
time history of the displacement on the ground, as seen in fig. 1201 In fact, we see in this figure 
that there are no substantial differences in the responses on the ground between the cases of the 
absence and presence of the block. 

8.2.5 Displacement in the substratum 

The focus here is on the displacement field in the subtratum when the solicitation is due to a 
deep source located at x* = (0,3000m). 

As shown in fig. 1211 the main components of the diffracted field (i.e. — n*) are those due to 
propagative and evanescent waves in the substratum, the latter being associated with a quasi-Love 
mode. 

The comparison of the results of the mode-matching and finite element methods is made in 
fig. 1221 The time history of the total displacement field is mainly composed of the incident and 
specularly-reflected flelds. Although a quasi-Love mode is excited, it hardly makes itself felt in the 
substratum. 
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Figure 20: Comparison of 27r times the spectrum (left panel) and of the time history (right panel) of 
the total displacement on the ground in the absence of the block (solid curves) and in the presence 
of the block (dashed curves) at various locations on the ground: (from the top to the bottom) 
at 150 and 300m from the center of a single 50m x 30m block. The shallow source is located at 
= (-3000,100m). 
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Figure 21: The solid curve represents 27r times the transfer function of the diffracted displacement 

at the interface between the substratum and the layer below the center of a single 50m x 30m block. 
The dotted curve represents the contribution of the interval Xi, the dashed curve the contribution 
of the interval Z2 and the dot-dashed curve the contribution of the interval J3. 




Figure 22: 2tt times the spectrum (left panel) and time history (right panel) of the total displace- 
ment at the interface between the substratum and the layer below the center of a single 50m x 30m 
block. The dashed curves result from the semi-analytical (mode-matching) computation, whereas 
the solid curves result from the numerical (finite element) computation . Note the splitting of the 
first peak. 
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9 Numerical results for two-block configurations in a Mexico City- 
like site 



Let us now consider the configuration involving two blocks, which is the simplest configuration 
for studying inter-block coupling effects. The latter go by the name of structure-soil-structure 
interaction. 

The two blocks, situated in a Mexico City- like environment, are located such that the center of 
the base segment of block 1 is (Om,Om) and that of block 2 is (-65m, Om). 

The parameters of the site are: = 2000 kg/m^, 0*^=600 m/s, = oo, = 1300 kg/m^, 
c^=60 m/s, = 30, with the soft layer thickness being /i = 50 m. The material constants of the 
blocks are: p"^ = 325 kg/m^, 0^=100 m/s, = iqo. 

The incident cylindrical wave is radiated by a deep line source located at (Om, 3000m). 

Recall that the eigenfrequencies of the displacement-free base block are UqJ^ = and 
the Haskell frequencies are v^^^^ = ^"^"^"^ wherein m = 0, 1, 2, .... Thus, the Haskell frequencies 
are 0.3, 0.9, 1.5 Hz 

If the zeroth-order quasi-mode coefficient is relevant, the dispersion relation of the configuration 
takes the from: 

(^(1) _ _ _^(2)^J _ _^(12)_^(21) ^^^^ (^33) 

wherein T)[ —T^ = 0, j = 1, 2 is the dispersion relation of the configuration with one block of char- 
acteristics of the block j (see eqEZJ and the term .7^2^^'' .7^2^^'' = ^cos {}?h\) Qqq^^^ (^cos [}?h'2) Qqo^^^ 
accounts for the coupling between the two blocks. 



9.1 Results relative to a two-block configuration consisting of a 50m x 30m block 
and a 40m x 40m block 

Block 1 is 50m high and 30m wide and block 2 is 40m high and 40m wide. Their center-to-center 
separation is 65m. 

The displacement-free base block eigenfrequencies are: v^^f^ = 0.5, 1.5Hz... for block 1, and 
z/^^s = 0.625, 1.875Hz... for block 2. 

Figure 1221 gives an indication of the frequency of occurrence of the modes of the system. An 
eigenfrequency is a frequency for which ?R. {T) = 0. The attenuation associated with a particular 
mode (at a frequency v*) is given by 

A quasi-Love mode is excited a,i v k, 0.3. We could expect two different quasi displacement- 
free base block modes for a system with two non-identical blocks (at i^qq^^ ~ 0.5 and 0.625-ff2;), 
but only one eigenfrequency, at ~ 0.7Hz, is found (actually this value is debatable, since the 
the eigenfrequency might actually correspond to a minimum of ||^|| rather than to a zero of J-). 
This means that: i) the modes of a complex configuration are not the union of the modes of the 
subsystems of which it is composed, ii) the single quasi displacement-free block mode results from 
the inter block coupling of the fields, i.e., the (structure- soil- structure interaction), iii) this mode 
is a coupled mode and is associated with a smaller attenuation than either of the corresponding 
modes of the two associated one-block configurations, iv) the correct expression of the modes of 
the system is the one given in section 11)31 involving the couplings between the blocks. In particular, 

(ij) 

the coupling matrix Qqq , i ^ j cannot be neglected. 
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Figure 23: Indications concerning the dispersion relation for a system of two blocks, one with 
dimensions 50m x 30m and the other with dimensions 40m x 40m. The solid curve represents 
and the dashed curve 

To distinguish this mode from the quasi displacement-free base block mode for a single block 
that accounts only for the structure-soil interaction, we call it the multi- displacement- free base block 
mode. 

9.1.1 Response on the top and bottom segments of the blocks 

Now consider the response at the centers of the top segments of the two blocks. The responses 
at these locations, as computed by the finite element and mode matching methods (with account 
taken only of the zeroth-order quasi-mode in each building), are seen in fig. [^to be almost identical 
in each block. The peak at VqI^^ ~ O.G-fTz, translates the excitation of the multi-displacement-free 
base block mode. This peak is sharper than the one encountered for only one block at z^qo^^ (due 
to excitation of the quasi displacement-free base block mode), which fact is essentially due to its 
larger amplitude. The important sharpness of this peak is also related to the fact that the excited 
mode is a coupled mode. This suggests that the larger the number of blocks, the larger will be the 
response at the resonance frequency of the multi-displacement-free base block mode. 

We compare in fig. [^the results obtained in (at the centers of the top and bottom segments of) 
the two blocks with those on the ground (at the same locations as the centers of bottom segments of 
the blocks) in the absence of these blocks. In the time domain, a small increase of the duration and 
a more substantial increase of the peak and cumulative amplitudes can be observed, in particular on 
the top segments of the two blocks. These increases are more important than in the configurations 
with a single 40m x 40m or 50m x 30m block (see figs. 1^ and ITHl . and is also somewhat more 
important in block 1 than in block 2, as already noticed in the cases of single blocks. The field 
vanishes at the center of the bottom segments at the displacement-free base frequencies: v = 0.5Hz 
for block 1 and = 0.625Hz for block 2. This emphasizes the fact that the multi displacement-free 
block mode is different from the quasi displacement-free block mode of the configuration with only 
one block, and is a manifestation of geometrical features. 
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Figure 24: 2it times the spectrum (left panels) and time histories (right panels) of the total displace- 
ment of a configuration with two blocks for a deep seismic line source located at x*' = (0,3000m). 
The top panels apply to the field at the center of the top segment of block 2 (40 x 40 block). 
The lower panels apply to the field at the center of the top segment of block 1 (50 x 30 block). 
The dashed curves result from a computation via the semi-analytical (quasi-mode) method (taking 
account only of the zeroth-order quasi-mode in each block), whereas the solid curves result from a 
finite element computation. 
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Figure 25: Comparison of 27r times the spectrum (left panels) and time history (right panels) of the 
total displacement on the ground in absence of blocks (solid curves) and in the presence of blocks 

(dashed curves). From the top to bottom: at the center of the summit segment of block 1, at the 
center of the summit segment of block 2, at the center of the base segment of block 1, and at the 
center of the base segment of block 2. The deep line source is located at x'^ = (0, —3000m). 
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Figure 26: Indications concerning the dispersion relation for a configuration with two 40m x 40m 
blocks. The sohd curve represents ?fi{T) and the dashed curve Q{J^)- 



9.2 Results relative to a two-block configuration with two 40m x 40m blocks 

The two blocks are 40m high and 40m wide and their center-to-center separation is 65m. 

The displacement-free base block eigenfrequencies of the two blocks are 0.625, 1.875 Hz ... 
Fig. 1261 gives an indication of the solutions of the dispersion relation = 0. The eigenfreqencies 
corresponding to ?R.J^ = 0, are seen to occur at 0.3, 0.35, 0.6, 0.7, 0.9, l.OHz... Some of them 
are close to each other and can be gathered together into the groups: 0.3 — 0.35, 0.6 — 0.7 and 
0.9 — 1.0 Hz. This suggests a splitting associated with the lifting of the degeneracy of eigenvalues. 
Nevertheless, the attenuations of the modes at 0.35, 0.7 and 0.9 — 1.0 are very large, so that we can 
expect them to be hardly excited (as will be confirmed in the following results). 

9.2.1 Response on the top and bottom segments of the blocks 

In fig. EZIwe see that the results obtained by the two computational methods pertaining to the 
response at the center of the summit segments of the two 40m x 40m blocks are nearly identical. 
The amplitude of the multi-displacement-free base block resonance peak is higher for two identical 
blocks, than it was for two different blocks. As pointed out during the previous discussion of the 
dispersion relation of the configuration, no splitting of the first and second peaks is noticed because 
of the high attenuation of one of each pair of split modes. On the contrary, the large response due 
to the excitation of the other one of the pairs of split modes is due to the small attenuation of these 
modes. Notice also that the second (quasi displacement-free base mode) peaks in the spectra of the 
figure are larger than the first peaks, in contrast to what was obtained for two dissimilar blocks. 

In fig. EHl we compare the displacement in the two blocks (in particular, at the centers of the 
summit and base segments) with the displacement on the ground in the absence of these blocks. 
The response of the two blocks are, in fact, identical, due to the deep source being located on the 
vertical dividing line between the two blocks. One observes an increase of the duration and of the 
peak and cumulative amplitudes that is larger than for two different blocks, in particular, on the 
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Figure 27: 27r times the spectra (left panel) and time histories (right panel) of the total displacement 

at the center of the summit segment of block 2 (top panels), and at the center of the summit segment 
of block 1 (bottom panels), as computed by the mode-matching method (accounting only for the 
zeroth-order quasi-mode) (dashed curves) and by the finite element method (solid curves) for a 
configuration with two 40m x 40m blocks solicited by the cylindrical wave radiated by a deep line 
source. 
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top segments. This increase seems to be due to the much more stronger response associated with 
the excitation of the multi displacement-free base mode block. This hypothesis is backed up by the 
fact that the said mode is more weakly excited on the base segment and at the same time the time 
history is much closer to that of the configuration with no blocks. Finally, the displacement at the 
center of the base segments vanishes at the frequency of occurrence of the displacement-free base 
mode of the block (which is not a mode of the global configuration) and is slightly different from 
the corresponding multi-displacement-free block mode eigenfrequency. 

9.3 Results relative to a two-block configuration with two 50m x 30m blocks 

The blocks are 50m high and 30m wide and their center-to-center separation is 65m. 

The displacement-free base block eigenfrequencies are 0.5, 1.5Hz, Fig. 1291 gives an indication 

of the solution of the dispersion relation. Eigenfreqencies (i.e., solutions of = 0) are found at 
0.3, 0.55, 1.0, 1.4... The attenuations associated with the quasi-Love mode at 0.3Hz and the multi 
displacement-free base block mode at 0.55Hz are rather small. The eigenfrequencies are not close 
to each other, contrary to the case of two identical 40m x 40m blocks. 

9.3.1 Response on the top and bottom segments of the blocks 

The response at the centers of the summit segments of both blocks computed by the finite 
element method is compared in fig. EOl to the corresponding response computed by the mode 
matching method (with account taken only of the zeroth-order quasi-mode. Once again, we observe 
these responses to be almost identical. Moreover, the amplitude of the multi- displacement-free 
base block resonance peak is observed (in the same figure) to be higher for two identical blocks, 
than it was for two different blocks. 

In fig. I^we compare the displacements in the two blocks (i.e., at the centers of the summit and 
base segments) with the displacement on the ground (at the same points as occupied by the centers 
of the base segments of the blocks) in the absence of the blocks. The response of the two blocks 
are identical for the previously- mentioned reasons. A much larger increase of the duration and of 
the peak and cumulative amplitudes (particularly on the top segments) is obtained for the two 
identical blocks than for two different blocks. This increase seems to be due to the much stronger 
response at the frequency corresponding to the excitation of the multi displacement-free base block 
mode. Once again, the displacement at the center of the base segments vanishes at a frequency 
corresponding to the occurrence of the displacement-free base mode of the block. 

9.3.2 Response on the ground 

In order to get another grip on the phenomena that are produced when two blocks are present, 
we now focus on the displacement field at points on the ground outside of the blocks. 
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Figure 28: Comparison of 2tt times the spectrum (left panels) and of the time history (right panels) 
of the total displacement on the ground in the absence of blocks (solid curves) and in the presence 
of the blocks (dashed curves). From the top to the bottom: at the center of the summit segment 
of block 1, at the center of the summit segment of block 2, at the center of the base segment of 
block 1, and at the center of the base segment of block 2. Two 40m x 40m blocks solicited by the 
cylindrical wave radiated by a deep line source Igigated at = (0, —3000m). 




Figure 29: Indications concerning the dispersion relation = 0. The sohd curves represent 3f?(jr) 
and the dashed curves 




Figure 30: Itx times the spectra (left panel) and time histories (right panel) of the total displacement 
at the center of the summit segment of block 2 (top panels), and of the block 1 (bottom panels), as 
computed by the mode-matching method (with account taken only of the zeroth-order quasi-mode 
(dashed curves) and the finite element method (solid curves). 



57 



500 r 



§400 



;}300 
^200 



2it||u^<^'{x,v)| 
2it||u\x,v)|l 



100r 



u2<^'(x,t) 
u\x,t) 





•^-50 



0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 
v(Hz) 



-lOO.L 



20 40 60 80 100 120 
t(s) 



500 r 



§400 



:;:-3oo 

^200 
#100 

4 

500 r 



27t||u2<2)(x,v)|| 
27t||u\x,V)|| 



lOOr 





0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 

v(Hz) 



-lOO.L 



20 40 60, 80 100 120 

Us) 



§400 



1300 
?200 



271||U^<^'(X,V)| 

27t||u\x,v)|| 




lOOr 



I 50 



-50 



-lOO.L 



u^'^>(x,t) 
u^x.t) 



20 40 60, 80 100 120 
t(s) 



500 r 



§400 



;;^-300 

^200 
flOO 




2,:||u2(2)(x,v)|| 100- 
27i||u\x,v)|| _ 





u2<2)(x,t) 
uVx.t) 



-lOO.L 



0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 20 40 60 80 100 120 

V (Hz) f(s) 



Figure 31: Comparison of 27r times the spectrum (left panels) and of the time history (right panels) 
of the total displacement on the ground in the absence of blocks (solid curves) and in the presence 
of blocks (dashed curves). From the top to the bottom: at the center of the summit segment of 
block 1, at the center of the summit segment of block 2, at the center of the base segment of block 
1, and at the center of the base segment of block 2. Two 50m x 30m blocks solicited by the 
cylindrical wave radiated by a deep line source Igigated at = (0, —3000m). 




Figure 32: 2tx times the transfer function (solid curve) of the total displacement on the ground 
at Xg = (150m, Om) (left panel), and Xg = (300m, Om) (right panel) for a two 50m x 30m block 
configuration solicited by the cylindrical wave radiated by a deep line source located at x® = 
(Om, 3000m). The dotted curve represents the contribution along the interval Ti, the dashed curve 
the contribution along the interval T2, and the dotted-dashed curve the contribution along the 
interval T3. 

Fig. 1321 depicts 27r times the transfer function ti^(x, u;)/S'(a') resulting from a computation by 
the mode-matching method. The influence of two blocks on the contribution of the various types 
of waves traveling in the layer is close to the one we noticed when only one block is present (see 
fig. I16|) . The amplitude of the first peak is more important than that of the one in fig. 1161 

The comparison in fig. 1331 between the displacement in the presence of the blocks and in their 
absence shows a small increase of the amplitude and of the duration in the time domain, this 
depending non-linearly on the location on the ground. In the frequency domain, the amplitude of 
the first peak is larger than for a single block. 

10 Snapshots of the displacement fields for one- and two-block 
configurations in a Mexico City-like site 

In order to better visualize the excitation of the quasi-Love mode due to the presence of one 
or two blocks in a Mexico City-like site, we show, in figs. and ESI the snapshots at various 
instants, of the displacement field, for one 50m x 30m block, and two identical 50m x 30m 
blocks, respectively, in response to the cylindrical wave radiated by a deep line source located at 
X* = (Om, 3000m). We notice, in the layer regions of these figures, waves that are re-radiated from 
the base segments of the blocks and which evolve into a field with a series of nodes and anti-nodes 
characteristic of a sum of modes dominated by the quasi-Love modes. Coupling between the two 
blocks (structure-soil-structure iteraction) is also noticeable in fig. |2Slat t = 21s. 
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Figure 33: Comparison of 27r times the spectrum (left panels) and of the time history (right panels) 
of the total displacement on the ground at Xg = (150m, Om) (top panels), and Xg = (300m, Om) 
(bottom panels) in the absence of blocks (solid curves) and in the presence of the blocks (dashed 
curves). Case of two 50m x 30m blocks solicited by the cylindrical wave radiated by a deep line 
source located at x* = (0, —3000m) 
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t = 30s t = 36s t = 42s 




Figure 34: Snapshots at various instants, t = 12s, t = 18s, t = 24s, t = 30s, t = 36s, t = 42s, 
of the total displacement field for one 50m x 30m-block configuration, solicited by the cylindrical 
wave radiated by a deep line source located at x*=(0m,3000m). 
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t = 42s 




Figure 35: Snapshots at t = 21s and t = 42s of the total displacement field for two identical 
50m X 30m- block configuration, solicited by the cylindrical wave radiated by a deep line source 
located at x*=(0m,3000m). 

11 Conclusions and preview of the contents of the companion pa- 
per 

The response, to the cylindrical wave radiated by a line source located in the substratum, of 
a finite set of non-equally sized, non-equally spaced blocks, each block modeling one, or a group 
of buildings, in welded contact with a soft layer overlying a hard half space, was investigated in a 
theoretical manner via the mode matching technique. 

The capacity of this technique to account for the complex phenomena provoked by the presence 
of blocks on the ground was demonstrated by comparison of the numerical results to which it leads 
to those obtained by a finite element method. 

It was shown that the presence of blocks induces a modification of the phenomena that are 
produced by the configuration without blocks, or of a configuration of closed blocks disconnected 
from the geophysical half-space. In particular, the blocks modify the dispersion relation of what, 
in the absence of the blocks, constitutes the Love modes. 

Three different and complementary points of view were developed (in the framework of the mode 
matching theory) concerning the dispersion relation relative to the configuration with blocks. The 
first emphasizes the role of the substructure (i.e. the geophysical (flat ground/layer /substratum) 
structure ). The second emphasizes the role of each particular block of the superstructure. The 
third point of view emphasizes the couplings of the fields in the superstructure with those in the 
substructure. 

The first two points of views enable the definition of two new types of modes relative to the 
configuration with blocks: the quasi-Love modes, which are small perturbations of the Love mode 
(which can exist when no blocks are present), and the quasi displacement- free base block modes, 
which are small perturbations of the displacement- free base block mode (which can exist when no 
geophysical structure is present). 
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These two types of quasi-modes account for coupling between a particular block of the super- 
structure with the substructure, but not for the couplings between blocks (when more than one 
block is present). Thus, the so-called soil- structure interaction is shown to be due to the excitation 
of quasi Love modes. 

The third point of view emphasizes the coupling between blocks (when more than one block 
is present). It was shown that this coupling, manifested by the existence of coupling matrices in 
the expressions for the quasi mode amplitudes, is of the same form as the coupling between one 
particular block and the substructure, which underlines the fact that the coupling between blocks 
is carried out via the substructure. 

The study of the dispersion relation for the multi-block configurations is very complex, but 
can be carried numerically. This was done for a two-block configuration and showed that a multi 
displacement-free base block mode can be produced. The latter was shown to correspond to a 
coupled mode, constituted by a combination of quasi displacement-free base block modes of each 
block, which are no longer excited as such in a configuration with more than one block. The multi 
displacement-free base block mode was shown to account for the so-called structure- soil- structure 
interaction. 

It was underlined that the modes of a complete one- or multi-block configuration do not con- 
stitute the reunion of the modes of the individual component structures. Thus, the phenomena for 
a complete N block configuration are not the sum of the phenomena for N — 1, N — 2, ..,0 block 
configurations. 

The excitation of these modes was then studied in the particular case of one and two blocks. 
A common feature of the infiuence of one/or more blocks is the excitation of the quasi-Love mode, 
which occurs even for solicitation by the waves radiated by a deep source (recall that, for this type 
of solicitation, it is not possible to excite ordinary Love modes in a fiat ground (i.e., no blocks) /soft 
horizontal layer /hard substratum configuration !23). The trace of quasi-Love mode excitation in 
the frequency domain was shown to be: i) for a deep source, a shift to lower frequency and an 
increase of the amplitude of the first (lowest-frequency) peak of the response, and ii) for a shallow 
source (the case in which waves whose structure is close to that of Love waves already exist in the 
layer and substratum in the absence of the block(s)), an increase of the amplitude, and little or no 
shift, of the first resonance peak. 

As concerns deep sources, the change of the phenomena, from a configuration without blocks 
(for which there exist only bulk waves in the geophysical structure), to one with blocks (for which 
there exist quasi-Love modes characterized by a field in the substratum that is predominantly a 
surface wave in the substratum) is a manifestation of the so-called soil- structure interaction. 

Quasi displacement-free base block modes are encountered only for configurations with one 
block whereas multi displacement-free base block modes are excited in configurations with two 
(and presumably, more) blocks. Both of these modes are excited with incident waves that are 
radiated either from deep or shallow sources. The multi displacement-free base block mode is 
more energetic and corresponds to a coupled mode. Its excitation is tied up with the so-called 
structure- soil- structure interaction. 

The modifications of frequency-domain response (and less so for the time-domain response) 
were found to be fairly substantial, even when only one block is present, at large distances from 
the block on the ground. 

Generally speaking, the peak and cumulative amplitudes and duration of the time histories are 
larger for a one- and two-block configuration than for flat ground, especially at locations within 
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the block, but no spectacular effects, such as those noticed during the Michoacan earthquake were 
found in the numerical simulations, either for the one- or two-block "cities". 

A configuration involving a larger number of blocks is difficult to study theoretically in the 
framework of the mode-matching formulation, which is why configurations with an infinite number 
of blocks, each one of which is close to the average shape and composition of the blocks of typical 
cities, is investigated in the companion paper. It is also shown in this paper that, beyond of the 
order of ten identical blocks, the response of the laterally-finite city is very close to that of the 
city with an infinite number of blocks. More importantly, it is found in the companion paper, that 
the phenomena are quite different, for configurations with a large number of identical blocks, from 
those for a small number of identical or non-identical blocks, and that the said phenomena for a 
large number of blocks are evocative of those observed during earthquakes in urban areas such as 
Mexico City. 

This indicates that it is probably not sufficient to carry out experiments and numerical simu- 
lations on an isolated building 1^ in order to predict correctly the seismic response of the 
building when it is surrounded by other buildings (the situation of most buildings in a typical city) . 

A An auxiliary problem: fields in the layer and in the substratum 
when the seismic distrubance takes the form of a ribbon source 
of width w located in the layer 

The problem is to determine the total fields n^(x,a;) and u^{x,uj), when the configuration 
without blocks (i.e. involving only a soft layer and a hard half space) is solicited by a horizontal 
ribbon source of width w centered at x** in the soft layer. 

A.l Boundary conditions and field representations 

The fields satisfy the continuity of displacement and traction conditions on Th, i.e., (|2l-ij) and 
H24|) . the vanishing traction condition on Tg H19|) . and the outgoing wave (radiation) condition. 
In the layer and in the substratum the field representations are: 

(^i(A:i,cu)e-i'=2^2 +5i(A:i,cu)e^'=2^M e^'^i^i^ , (139) 

and 

u°(x,u;) = r i30(A;i,u;)ei(^i-i+'=§-2)^ ^ (140) 

J— oo ^2 

respectively. 

A. 2 The (incident) field radiated by a ribbon source of width w 

By employing the Green's theorem, the incident field can be written as: 

u'\^,oj)= [ G(||x-y||,c^)s^i(y,c^)dzi7(y); Vx,eM^ (141) 

JR2 
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wherein y is a vector from the origin O and pointing to a generic point in the sagittal plane (yi, 2/2) 
G(||x — y||,a;) the 2D free-space Green's function, and s*-^(y, cj) the source density function: 



/i(y,a;) = ^[F(yi-x? + |)-/?(yi-x?-|) 8 {y2 - x?,) , (142) 



w 



wherein H{ . ) is the Heaviside function and 5{ . ) the Dirac delta distribution. 
Introducing (|142|) into (|141j) . enables the incident field to be written as: 



u 



*i(x,a;) = i- rSMei('=i(^^-^i)+'=Jl^^-^^l)sincfA:i-') ^ , (143) 
4vr7oo \ 2J kh 



wherein sinc((^) := 

A. 3 Expression of the fields in the the presence of a ribbon source of width w 

The introduction of the field representations (|139|) and (|14()|) into the boundary conditions (with 
the appropriate projection), leads, after the resolution of the resulting linear system, to: 

• for X2^X2'- 

1 , , \S{uo)^i^ r "^^^ {A - h)) + i|S sin {kl (xi - h)) 

u [x,uj) = — / —J 

2vr/iU J-00 cos (klh) -i^ sin {klh) 

cos {klx2) ei^^(--^1)sinc (k,^) ^ , (144) 

K2 



for X2 S 



1 , , iS(u;)f,^ r "^^^ (^2 i.X2 - h)) + i|^ sin {kl [x2 - h)) 
ti (x, w) — ' 



2vr^" 7-00 cos(A:i/i) -i^sin(A;i/i) 

w\ dki 



cos (fcixl) e^^'^^'-^ihinc ^ , (145) 



and 



u^'ix.uj) = — / -1 smc fci— cosffcoa^oj-TTr • (146) 

^ ' 2vr^° J -00 cos (klh) - i^i sin (klh) ^ 2 ^ ^ ^ 2; , ^ ^ 



-00 cos (A;2/i) - i^^fjp sin (A;2/i) V 2/ ^2 
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